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FOREWORD 

This  report  was  prepared  for  the  United  States  Air  Force  by  the 
Cornell  Aeronautical  Laboratory,  Inc.,  Buffalo,  New  York  as  required  by 
Contract  AF33{&16)-?572.  v  . 

Nl 

The  work  reported  herein  was  performed  by  the  Flight  Research  De¬ 
partment  of  Cornell  Aeronautical  Laboratory  under  the  sponsorship  of  the 
Flight  Control  Laboratory  of  the  Aeronautical  Systems  Division,  Wright  - 
Patterson  Air  Force  Base,  Dayton,  Ohio,  as  Project  No.  8225;  Task  No. 
82181,  to  perform  research  into  basic  considerations  of  adaptive  flight  con¬ 
trol  systems.  The  projectwas  administered  by  Lt.P.C.  Gregory  of  the  Flight 
Control  Laboratory. 
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on  this  project  were  R.G.  Buscher  of  the  General  Electric  Company,  D.O. 
Dommasch  of  DODCO,  Inc. ,  and  G.  W.  Anderson  of  Aeronutronic. 
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ABSTRACT 

This  report  contains,  in  two  volumes,  the  results  of  a  study  of  the 
application  and  evaluation  of  certain  adaptive  control  technique ^1^8  applied  to 
advanced  vehicles  of  the  X-15  and  Dyna-Soar  type^  Volume  I  deals  with  those 
^^aspects  of  the  study  pertaining  to  the  General  Electric  Self-Adaptive  Flight 
^  Control  System,  Volume  II  deals  with  the  DODCO,  Inc.,  Optimum  Limited- 
V  Information  Adaptive  Flight  Control  System. 


Voliune  I^The  general  concept  of  adaptive  control  through  the  use  of 
reference  models  is  discussed,  and  particular  models  are  evaluated  based  on 
the  present  status  of  airplane  handling  qualities  research.  The  G.E.  System 
is  applied  to  the  problem  of  controlling  the  longitudinal  short-period  motions 
of  the  X-15  airplane  during  re-entry.  -~The  first  treatment  of  the  problem  is 
rather  general,  and  uses  essentially  linear  techniques  to  investigate  die  ref¬ 
erence  model  concepts,  selection  of  system  parameters,  responses  to  command 
inputs  and  gusts,  effect  of  basic  airplane  static  and  dynamic  instability,  effect 
of  sensor  dynamics,  and  dynamics  of  the  adaptive  loop.  Then  a  more  detailed 
study  is  made  of  certain  problem  areas,  ^including  the  response  to  inputs  at 
die  actuator  valve,  effect  of  actuator  non-linearities,  frequency  sensor  charac¬ 
teristics,  and  the  effect  of  noise  in  the  primary  sensor,  in  this  case  a  rate  gyro.««^^_^ 
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heading  angle,  angle  between  reference  azimuth  (North) 
and  the  projection  of  the  x  axis  in  the  horizontal  plane 

©  attitude  angle,  angle  between  x  axis  and  the  horizontal  plane 

0  bank  angle,  angle  between  y  axis  and  a  horizontal  line  in 
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Control  System  and  Transfer  Function  Notation 

A  transfer  function  is  defined  for  an  element  or  system  as  the  ratio 
of  the  Laplace  transforms  of  output  to  input.  For  a  linear  system,  this 
represents  the  transform  of  the  impulse  response. 

t  time 

•  d(  ) 

(  )  time  derivative,  - 

3  complex  frequency  S»a+j®,  Laplace  transform  variable 

«  . 

F(s)  Laplace  transform  of  f(t),  F(s)  *  f  meat 

G( s)  forward  loop  transfer  function 
Hfc)  feedback  loop  transfer  function 
P(s )  prefilter  transfer  function 
N(s)  numerator  of  transfer  function 
Dfs)  denominator  of  transfer  function 
Mfs)  reference  model  transfer  function 
C  command  input 

$  response  output 

D  disturbance  input 
6  error  signal 
6)  frequency,  rad/sec 
£*  damping  ratio 
6)4  undamped  natural  frequency 
time  constant 
f  frequency.,  cps 
7"  period 
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SECTION  I 
INTRODUCTION 

1.  1  Philosophy 

The  automatic  control  specialist  faces  a  serious  problem  when  trying 
to  design  a  flight  control  system  for  a  modern,  high-speed  flight  vehicle.  Vehicle 
characteristics  vary  so  much  that  it  is  extremely  difficult  to  design  a  satisfactory 
constant  parameter  control  system.  The  concept  of  using  feedback  in  a  control 
system  was  originally  devised  to  overcome  such  changes  in  system  characteristics. 
However,  for  supersonic  and  hypersonic  vehicles,  these  variations  are  so 
extreme  that  the  control  system  designer  finds  it  impossible  to  devise  satisfac¬ 
tory  constant  parameter  control  systems.  He  is  forced  to  program  the  control 
system  parameters  as  functions  of  the  flight  condition  —  for  example,  a  feedback 
gain  as  a  function  of  dynamic  pressure.  To  successfully  design  such  a  control 
system,  the  designer  usually  requires  extensive  knowledge  of  the  vehicle's 
characteristics,  and  often  sufficient  data  is  not  available  until  after  the  vehicle 
is  flight  tested. 

Recently,  considerable  attention  has  been  focused  on  the  "self-adaptive" 
control  concept.  Self-adaption  implies  measurement  of  the  system's  dynamic 
response,  evaluation  of  this  response  by  some  criteria  or  by  comparison  with  a 
desired  response,  and  adjustment  of  the  parameters  of  the  control  system  to 
bring  the  measured  response  into  agreement  with  the  desired  response.  This 
is  essentially  the  process  used  by  the  control  system  designer,  but  a  self- 
adaptive  system  performs  the  process  itself.  Basically,  an  adaptive  system 
is  one  in  which  the  control  system  parameters  are  automatically  adjusted  to 
compensate  for  changing  vehicle  characteristics,  either  due  to  a  changing 
environment  or  due  to  internal  changes;  a  self-adaptive  system  is  one  in  which 
this  adjustment  is  made  by  internal  measurement  and  evaluation  of  the  dynamic 
response  itself  rather  than  by  measurement  of  some  environmental  characteristic. 
An  environmental  characteristic  is  one  which  is  only  indirectly  correlated  with 
the  response  of  the  control  system  being  considered.  An  internal  process  is 

Manuscript  released  by  author  31  July  1961  for  publication  as  an  ASD  Technical 
Report. 
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one  which  is  included  in  the  control  system  directly,  and  internal  changes  are 
correlated  with  the  response  of  the  control  system  being  considered.  The 
terms  "adaptive"  and  "self- adaptive •'  are  conceptual  in  nature  as  is  the  term 
"feedback",  and  their  implication  is  not  on  what  quantity  is  measured  but  rather 
why  it  is  measured  and  what  is  done  with  the  measurement. 

To  see  what  adaptive  control  may  have  to  offer,  the  problem  must  be 
formulated  more  precisely.  The  object  is  control  the  dynamic  performance  of 
a  flight  vehicle  in  some  desired  maimer.  The  vehicle's  characteristics  change, 
and  the  change  is  so  extreme  that  conventional  feedback  control  techniques  will 
not  suffice  to  compensate  adequately  for  this  change.  Accordingly,  compensating 
changes  could  be  programmed  in  the  control  system  parameters  as  a  function  of 
the  environmental  condition  that  produces  the  change.  The  real  difficulty  lies 
in  obtaining  a  sufficiently  accurate  description  of  the  vehicle's  characteristics 
as  a  function  of  the  environment.  What  is  needed  are  methods  for  simplifying 
the  design  problem  by  minimizing  its  dependence  on  information  concerning  the 
vehicle' 8  characteristics.  There  are  a  number  of  possibilities.  Part  of  the 
design  task  can  be  shifted  to  the  control  system.  New  techniques  can  be  developed 
to  measure  the  dynamic  performance  directly,  and  use  this  information  to  adjust 
the  control  system  parameters  without  recourse  to  knowledge  of  the  vehicle's 
characteristics.  New  types  of  control  system  elements  —  presumably  nonlinear 
ones  —  can  be  developed  which  would  be  insensitive  to  changes  in  the  vehicle's 
characteristics.  More  exotic  concepts,  such  as  thinking  or  self-organizing 
systems,  implying  drastic  changes  in  control  philosophy,  can  be  readily  con¬ 
ceived  though  their  implementation  represents  a  formidable  task.  However, 
practically  speaking,  it  can  be  assumed  that  the  basic  tenets  of  feedback  control 
theory  will  be  retained,  and  that  improvements  in  the  design  of  flight  control 
systems  can  be  expected  to  result  from  the  introduction  of  logical  or  decision¬ 
making  elements,  non-linear  elements,  direct  measurements  of  dynamic  per¬ 
formance,  and  other  similar  innovations.  Clearly,  the  introduction  of  new 
techniques  and  the  concepts  of  self-adaptive  control  offer  methods  for  designing 
better  flight  control  systems. 
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Reflection  on  the  principal  purpose  for  introducing  adaptive  control 
(and  also  simple  feedback  control)  leads  to  a  rather  interesting  and  significant 
conclusion.  The  purpose  is  to  provide  for  a  wider  latitude  of  vehicle  or  plant 
characteristics,  or,  to  minimize  the  need  for  knowledge  of  these  characteristics. 
Consider  two  extreme  hypothetical  cases:  one,  the  case  where  the  plant  char¬ 
acteristics  are  known  completely;  two,  the  case  where  only  the  input  and  output 
quantities  are  known.  Further,  consider  that  a  particular  and  explicitly  defined 
response  is  desired,  and  that  the  sole  basis  for  performance  evaluation  is  the 
error  between  the  achieved  response  and  the  desired  response.  To  obtain  the 
desired  response  in  the  first  case  requires  only  an  open  loop  input  —  no  feedback 
and  no  adaption  —  and  in  theory  the  error  can  be  made  identically  zero.  To 
obtain  the  desired  response  in  the  second  case  requires  the  ultimate  in  adaption  - 
the  response  of  the  plant  to  the  input  must  be  measured,  interpreted,  and  then 
used  to  modify  the  input  (compensation),  and  this  must  be  done  continually. 
Because  knowledge  of  the  system  can  only  be  achieved  after  some  finite  response 
has  been  measured,  and  because  the  measurement,  evaluation,  and  compensation 
process  requires  a  finite  time,  there  must  necessarily  always  be  a  finite  error. 
Thus  in  theory,  the  first  system  will  always  outperform  the  second  system, 
because  it  is  designed  with  complete  a  priori  knowledge  of  the  plant  while  the 
second  system  must  make  an  identification  of  these  dynamics.  Generalizing 
leads  to  a  rather  significant  conclusion.  The  use  of  an  adaptive  process  leads 
to  a  degradation  or  slowing  up  of  the  control  process;  however,  performance 
gains  can  be  achieved  through  the  indirect  effects  of  adaption,  that  is,  by  making 
the  system  less  sensitive  to  unknown  or  changing  characteristics.  This 
phenomenon  is  commonly  demonstrated  in  digital  computer  programs,  where  it 
is  well  known  that  the  addition  of  logical  or  decision  steps  slows  the  computing 
process.  Thus,  just  because  a  system  is  adaptive  does  not  mean  it  will  perform 
better,  and  making  a  system  adaptive  where  it  is  not  required  may  be  actually 
undesirable.  Full  advantage  should  be  taken  of  all  available  knowledge  of  the 
components  in  designing  any  control  system. 
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1.2  Purpose 

A  study  of  some  basic  considerations  of  adaptive  flight  control  systems 
has  been  undertaken  by  the  Flight  Research  Department  of  the  Cornell  Aeronau¬ 
tical  Laboratory  for  die  Flight  Control  Laboratory  of  the  Aeronautical  Systems 
Division,  under  the  direction  of  Lt.  Philip  Gregory.  Some  of  the  earliest  work 
in  this  field  was  performed  by  Campbell  of  the  Cornell  Aeronautical  Laboratory 
(References  1  and  2),  and  interest  in  this  work  has  continued  since.  The  Flight 
Control  Laboratory  has  sponsored  considerable  research  on  adaptive  flight  con¬ 
trol  systems,  and  gave  impetus  to  much  of  the  earlier  work  (Reference  3).  The 
results  have  led  to  the  formulation  of  a  number  of  different  schemes  for  adaption. 
During  the  course  of  the  present  project,  three  of  these  schemes  -  ones  devel¬ 
oped  by  General  Electric,  DODCO,  and  Aeronutronic  -  were  evaluated  with 
respect  to  their  application  to  advanced  flight  vehicles  of  the  X-15  and  Dyna-Soar 
type.  The  intent  of  this  study  was  not  to  determine  which  system  would  be  best  for 
the  X-15  or  Dyna-Soar;  rather,  the  intent  was  to  achieve  a  more  thorough  under¬ 
standing  of  self-adaptive  control  systems  -  how  do  they  behave?  how  can  they  be 
designed?  what  problems  will  be  encountered?  what  kind  of  adaptive  elements 
are  most  likely  to  solve  these  problems  ?  The  study  was  directed  toward  the  prob¬ 
lems  associated  with  hypersonic  flight  vehicles,  but  the  results  are  generally 
applicable  to  the  control  of  any  flight  vehicle  whose  characteristics  vary  over  a 
wide  range. 

1 . 3  Approach 

The  approach  used  in  this  study  may  be  outlined  briefly  as  follows.  First, 
a  detailed  description  and  specification  of  the  airframes  and  the  adaptive  systems 
was  needed.  Obtaining  the  airframe  data  was  no  problem  since  data  for  the  X-15 
and  Dyna-Soar  were  already  available  in  house.  But  for  the  adaptive  control 
systems,  not  only  was  it  necessary  to  obtain  all  the  available  published  material, 
but  it  was  also  necessary,  for  a  number  of  reasons,  to  contact  the  people  re¬ 
sponsible  for  the  development  work  on  each  of  the  three  systems,  and  to 
establish  a  working  relationship  with  these  people.  None  of  the  three  systems 
was  sufficiently  well  defined  in  the  published  literature,  because  the  development 
work  was  very  recent  or  still  underway  and  because  of  the  proprietary  nature  of 
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some  of  this  work.  Furthermore,  it  was  desirable  not  only  to  obtain  an  accurate 
physical  description  of  each  of  the  three  systems,  but  also  to  obtain  a  good  under¬ 
standing  of  the  underlying  design  philosophy.  If  the  evaluation  of  the  capabilities 
of  the  three  systems  was  to  be  meaningful,  each  would  need  the  full  benefit  of  the 
developer's  work  and  the  best  thoughts  and  concepts  concerning  their  application. 

A  general  plan  was  used  to  study  the  various  systems.  First,  the  adaptive 
concepts  were  applied  under  idealized  conditions  -  optimum  conditions  —  to  obtain 
a  measure  of  the  best  performance  that  could  be  realized.  Practical  considerations 
were  next  introduced  to  see  what  degradation  in  performance  occurred.  Finally, 
the  over-all  scope  of  the  work  was  broadened  to  include  n?w  considerations.  It 
was  not  intended  to  expend  equal  effort  on  each  of  the  three  systems  or  the  two 
vehicles  being  considered.  The  effort  was  expended  where  it  could  be  expected 
to  yield  the  most  results,  and  the  direction  of  this  effort  was  re-oriented  when¬ 
ever  it  was  deemed  appropriate. 

The  airframe  characteristics  specifically  studied  were  restricted  to 
longitudinal  motions.  Speed  was  not  included  as  a  degree  of  freedom,  but  as  an 
independent  variable  or  a  parameter.  The  equations  of  motion  for  the  airframe 
were  linearized.  In  application  of  the  adaptive  control  systems,  only  pitch-rate 
signals  were  used  for  control,  though  other  variables  (flight  path  angle,  angle 
of  attack,  etc. )  were  monitored.  Sensor  and  actuator  dynamics  were  included. 

Both  linear  and  non-linear  analysis  techniques  were  used.  The  root  locus 
technique  was  relied  upon  heavily  for  the  linear  analysis  because  an  ESIAC 
computer  was  available  (described  in  Appendix  B).  Analog  computer*  were 
used  for  computing  actual  responses  —  for  convenience  in  the  linear  case,  of 
necessity  in  the  non-linear  one. 

When  the  required  data  describing  the  three  systems  and  the  two  air¬ 
frames  was  in  hand,  a  preliminary  evaluation  was  made  of  the  problems  and 
difficulties  associated  with  applying  each  of  the  systems  to  the  two  vehicles. 

This  evaluation  led  to  the  following  conclusions.  The  X-15  would  be  suitable 
for  most  of  the  work,  since  it  represented  a  defined  configuration,  and  em¬ 
bodies  most  of  the  longitudinal  short  period  problems  that  would  be  encountered 
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in  Dyna-Soar.  The  G.E.  system  was  the  best  defined  and  its  principle  of  oper¬ 
ation  the  most  easily  understood  of  the  three  systems,  and  hence  would  be  the 
first  analyzed.  The  DODCO  system,  though  not  so  well  defined  or  as  easily 
analyzed  as  the  G.E.  system,  was  still  a  complete  concept  and  would  be  suit¬ 
able  as  the  second  candidate  for  study.  The  Aeronutronic  system  was  not  a 
complete  concept.  Rather,  it  represented  a  developed  except  for  measure¬ 
ment  of  dynamic  performance;  the  application  of  the  measurement  concept  to 
adaptive  control  still  required  further  research.  In  addition,  actual  hardware, 
in  the  form  of  a  digital  cross-correlator,  would  have  to  be  used  to  obtain  ana¬ 
log  simulation.  Thus,  study  of  the  Aeronutronic  concept  would  best  be  deferred. 
A  brief  description  of  the  concept  of  each  system  and  how  it  was  evaluated  is 
given  ir.  the  succeeding  paragraphs. 

The  G.E.  system  (Reference  4)  embodies  a  frequency  sensor  as  the 
basic  adaptive  element  and  a  reference  model  in  the  feedback.  Conceptually, 
the  frequency  sensor  maximizes  loop  gain  while  the  reference  model  holds  the 
response  constant.  These  concepts  were  applied  to  synthesize  a  pitch- rate 
control  system  for  the  X-15,  and  root  locus  techniques  and  analog  computers 
became  the  primary  tools  for  system  ana’ysis  and  evaluation.  Among  subjects 
studied  were  the  reference  model  concepts  sensor  dynamics,  actuator  dynamics 
and  non-linearity,  adaptive  loop  dynamics,  and  frequency  sensor  characteristics. 

The  DODCO  System  (References  5,  6,  and  24)  employs  in  the  main  con¬ 
trol  loop  a  dual-gain  adaptive  element  -  a  gain  switching  device  -  which  is 
controlled  by  functions  of  the  error  and  error  rate.  Here  the  error  is  defined 
as  the  difference  between  the  actual  response  as  measured  and  the  desired  re¬ 
sponse  as  obtained  from  a  reference  model.  The  system  is  non-linear;  though 
simple  in  mechanization,  it  presents  the  analyst  with  a  formidable  problem. 

The  key  concept  lies  in  the  switching  logic  which  derives  from  asking  the  ques¬ 
tion:  What  is  the  optimum  control  policy  for  ininimizing  the  error?  The  sys¬ 
tem  as  applied  to  pitch-rate  control  was  mechanized  on  an  analog  computer, 
and  the  first  task  was  to  duplicate  some  of  the  results  in  Reference  5  which  were 
obtained  using  a  digital  computer.  Then  system  parameters  were  varied  on  the 
analog  computer  to  evaluate  performance  and  stability.  Emphasis  was  placed 
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on  determining  the  limitations  of  the  system  and  on  developing  synthesis 
methods  for  use  in  design. 

The  Aeronutronic  system  (Reference  7)  embodies  the  principle  of  using 
die  impulse  response  of  the  vehicle,  obtained  from  the  cross -correlation  of 
the  output  to  a  tailored  white  noise  input,  to  adapt  the  flight  control  system. 
Study  of  the  control  concepts  led  to  the  following  conclusions.  Though  Aero- 
nutronic  had  demonstrated  a  technique  for  generating  the  impulse  response  and 
had  conceived  a  method  for  controlling  the  damping  ratio  of  a  second  order  sys¬ 
tem,  no  satisfactory  way  of  adapting  a  flight  control  system  using  the  impulse 
response  was  readily  available.  It  appeared  that  development  of  such  an  adap¬ 
tive  technique  would  be  feasible  and  practical,  but  would  require  substantial 
effort.  One  limitation  appeared  to  be  clearly  present:  there  is  an  important  and 
fundamental  time  delay  involved  in  the  adaptive  process;  this  time  delay  is  the 
time  it  takes  to  obtain  the  impulse  response,  and  must  be  at  least  as  long  as  the 
time  it  takes  for  the  impulse  response  to  decay  to  sensibly  zero.  In  view  of  the 
difficulties  involved  and  the  limited  effort  available,  further  study  of  the  Aero- 
nutronic  concept  was  abandoned. 
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SECTION  II 

DEFINING  THE  CONTROL  PROBLEM 

Z.  1  General  Considerations 

The  control  problem  for  a  hypersonic  or  re-entry  vehicle  is  a  multi- 
variable  one,  involving  many  modes  of  motion.  The  pilot  or  flight  director 
must  control  the  vehicle's  attitude  with  respect  to  the  flight  path  (at  and ), 
its  orientation  in  space  ( ,  &  ,  and  ),  and  its  flight  path  (  h  ,  /  ,  V 

and  flight  path  azimuth).  Control  is  normally  obtained  through  four  vehicle 
control  inputs:  the  pitch,  roll, and  yaw  controls  either  of  the  aerodynamic  or 
reaction  type,  and  the  thrust.  It  is  quite  obvious  that  independent  control  of 
all  the  output  variables  is  not  possible.  Usually,  attitude  control  is  associated 
with  the  rapid  modes  of  motion;  flight  path  dynamics,  with  the  slow  modes  of 
motion.  Accordingly,  the  control  problem  and  control  systems  are  broken  up 
into  two  general  groups,  which,  following  Dommasch  (Reference  5),  may  be 
labeled  as: 

1.  The  micro-control  problem  -  associated  with  the  fast  modes 
of  motion  and  what  are  often  called  the  "inner"  loops  in  the 
control  system  -  and 

Z.  The  macro-control  problem  -  associated  with  the  slow  modes 
of  motion  and  what  are  often  called  the  "outer"  loops  in  the 
control  problem. 

The  micro-control  problem  has  two  objectives:  (1)  to  provide  Satisfactory 
dynamics  for  the  far.t  modes  of  motion,  and  (2)  to  provide  suitable  responses 
of  the  fast  motions  for  use  in  the  macro-control  problem.  The  second  objective 
implies  providing  the  proper  statics  and  insuring  that  the  fast  modes  of  the 
micro-control  problem  are  fast  enough  to  avoid  interaction  with  the  slow  modes 
of  the  macro-control  problem.  Of  course,  such  interaction  cannot  always  be 
avoided,  but  it  can  be  minimized. 
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In  the  present  study,  we  are  concerned  with  only  the  micro-control 
problem  and  with  only  the  longitudinal  aspects  of  that  problem.  Thus  the  concern 
is  to  provide  a  proper  longitudinal  short-period  response  for  the  airplane  — 
the  micro-control  problem  -  and  insuring  that  the  adaptive  process  is  fast  enough 
to  account  for  the  changes  associated  with  flight  condition  -  the  macro-control 
problem.  More  severe  requirements  on  the  adaptive  process  may  be  posed  by 
changes  that  occur  rapidly  -  perhaps  with  angle  ot  attack  or  due  to  some  control 
system  non-linearity  such  as  backlash. 

The  general  approach  to  the  control  problem  may  be  listed  in  several 

steps: 

1.  Defining  the  vehicle's  characteristics, 

2.  Defining  a  desired  response, 

3.  Developing  general  control  schemes  which  will 
produce  the  desired  response, 

4.  Determining  the  problems  associated  with  the  control 
schemes  and  the  adverse  conditions  that  they  will  be 
subjected  to, 

5.  Implementing  the  control  schemes  into  concrete  systems,  and 

6.  Evaluating  the  efficacy  of  these  systems  under  expected 
adverse  conditions. 

This  study  is  primarily  concerned  with  the  last  step,  but  from  a  theoretical 
standpoint  rather  than  a  hardware  standpoint.  However,  all  these  steps  are  in¬ 
volved,  though  the  vehicles  and  the  control  schemes,  and  even  their  implemen¬ 
tation  to  some  extent,  are  defined  and  are  the  starting  point  for  this  study. 

2.  2  High-Gain  Control  Schemes 

The  control  schemes  used  in  this  study  are  based  on  those  developed  in 
feedback  control  theory,  but  modified  to  provide  greater  adaption  to  changes  in 
the  elements  of  the  system. 
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The  general  single  loop  problem  -  one  input  and  one  output  -  can  be 
reduced  to  the  block  diagram  shown  in  Figure  1, 


FIGURE  1  BI<OCK  DIAGRAM  OF  SINGLE -INI- UT 
SINGLE-OUTPUT  FEEDBACK  SYSTEM 

C  -  command 
I  =  input  to  closed  loop 
R  =  response 

0  =  feedback 

£  =  I-£i  ,  the  error 

P  =  prefilter  transfer  function 

KG  =  forward  loop  transfer  function 
W  =  feedback  loop  transfer  function 

The  response  to  a  command  input  is 


_R_  -  P  /'Km  \ 
C  *  H ( I  +  K6H ) 

and 

L'ltn  R  _  P 
C  "  H 


(2.D 


(2.<S), 
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If  the  gain  K  can  be  maintained  at  a  high  value,  then  the  response  of  the 
system  will  be  approximated  by  P/  H  .  We  have  control  over  P  and  H  , 
and  hence  the  ability  to  tailor  the  response  of  the  system.  We  may  formalize 
this  with  the  following  equation: 

/_E\  .  l;m  £.  /KM  \ 

UJU,;*-*-  h  t <+*'««/  i2-3) 

The  difference  between  tailoring  the  response  with  a  prefilter  or  with 
feedback  control  has  two  basic  facets.  The  first  is  purely  a  practical  one. 

With  prefiltering  only,  and  unity  feedback 

-  P  (2-4> 

v  ^  'DESIRED 

With  feedback  compensation  only, 

and  H  is  the  inverse  of  the  desired  response.  In  practice,  it  may  be  difficult 
to  generate  one  of  the  two  forms  -  most  likely  the  inverse.  The  second  differ¬ 
ence  between  the  two  methods  of  compensation  is  a  more  fundamental  one.  Con¬ 
sider  the  case  shown  in  Figure  2  ,  where  a  disturbance  enters  the  system  at 
some  interior  point  in  the  forward  loop  transfer  function. 


D 


FIGURE  2  BLOCK  DIAGRAM  OF  DUAL-INPUT 
SINGLE-OUTPUT  FEEDBACK  SYSTEM 
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K&  -  K;  Q,  Kz  Gz 

P  =  Disturbance  Input 

The  response  to  command  inputs  is  unchanged  (Equation  2.  1),  but  the  response 
to  a  disturbance  is  given  as 


JL  /  /  KCH  \ 

D  c*o~  + 


For  high  loop  gain, 

JL  £  __J _ 

D  K,&,H 


(2.6) 


(2.7) 


If  KL  is  finite  and  the  high  loop  gain  is  obtained  by  making  Kf  large,  then 


Li  hn 


(2.8) 


Thus,  the  important  difference  between  prefilter  type  compensation  and  feedback 
compensation  lies  in  their  different  effect  on  the  response  to  disturbance  inputs. 
These  disturbance  inputs  may  be  external  as  from  gust  inputs,  or  internal  as 
from  noise  in  electronic  components;  or  they  may  represent  cross-coupling  inputs 
from  other  portions  of  the  control  system.  In  general,  one  attempts  to  minimize 
the  response  to  disturbances.  Thus,  it  may  be  advisable  to  tailor  R/0  with  H  » 
and  ^  ,  while  P  is  used  to  tailor  R/C  . 

fX\  *  JL 

Desired  H  (2'9) 

/JL\  -  <  <2'10> 

^  D  ^OESiRBD  K  G,  H 
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j  The  role  of  Kg  Gg  is  of  significance,  and  can  be  illustrated  by  rewriting 

A  Equation  2.6  as 


*£* 


h 


-Hi  =jk^_ 

D'c.0  l  +  KGH 


(2.11) 


The  disturbance  response  is  proportional  to  K^  .  If  is  unity,  then  the 

disturbance  effectively  occurs  directly  in  the  output;  and  the  higher  the  fre¬ 
quency  content  of  the  disturbance  (beyond ‘a  certain  frequency),  the  less  will  be 
the  effect  of  the  control  system  in  minimizing  the  disturbance.  For  instance,  a 
step  in  D  will  result  initially  in  a  Mtep  in  R  . 


X 
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So  far  the  concept  of  adaption,  in  the  sense  of  adjustment,  has  not  been 
introduced.  However,  adaption  in  the  broader  sense,  i.e.,  ability  to  cope  with 
changes  in  the  controlled  element  in  the  forward  loop,  has  been  considered. 

The  simple  feedback  systems  of  Figure  1  and  2  have  been  discussed  from  the 
adaptive  point  of  view.  That  is,  the  open  loop  gain  K  has  been  assumed  suf¬ 
ficiently  high  that  the  response  has  been  formulated  in  terms  of  P  and  H  , 
without  regard  to  the  specific  value  of  K  or  the  form  of  6  .  This  is  not  a  new 
approach  (e.  g. ,  Reference  9,  Chapter  6,  p.  322),  but  it  has  not  been  applied  in  flight 
control  systems  for  good  reasons.  The  value  of  the  loop  gain  K  contains  the 
elevator  effectiveness  which  varies  over  a  wide  range  (considering  longitudinal 
control  and  just  one  of  the  varying  parameters),  and  it  is  difficult  to  design  a 
control  system  with  a  sufficiently  high  constant  gain  to  accommodate  the  vari¬ 
ations.  The  dilemma  is  sluggish,  inadequate  response  at  one  end  of  the  spec¬ 
trum,  instability  at  the  other  end.  In  practice,  control  systems  have  been  de¬ 
signed  with  variable  gains  using  air  data  parameters  such  as  dynamic  pres¬ 
sure  for  adjustment.  The  results  have  often  been  unimpressive,  and  have  led 
to  the  quest  for  better  concepts  of  adaption. 


Applying  the  high-gain  adaptive  concept  to  feedback  control  in  an  optimum  manner 
implies  developing  techniques  for  keeping  the  ioop  gain  at  the  highest  value  com¬ 
patible  with  stability.  Most  of  the  porposed  adaptive  control  schemes  approach 
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the  problem  in  this  manner.  Some  measurement  of  the  response  is  made,  in¬ 
terpreted  in  terms  of  stability,  and  then  used  to  adjust  the  loop  gain.  Thus,  an 
organized  and  clear  understanding  of  the  linear  feedback  control  problem  is  a 
necessary  prerequisite  to  the  study  of  adaptive  control  systems. 

2.3  Model  Approach  to  Control  System  Design 

The  concept  of  using  a  reference  model  in  flight  control  systems  was 
probably  first  proposed  by  Campbell  (Reference  1).  The  concept  was  suggested 
by  studies,  reported  in  References  10,  11  and  12,  which  showed  that  an  "opti¬ 
mum"  response  for  the  longitudinal  short-period  mode  could  be  specified  in 
terms  of  a  second-order  model  with  constant  characteristics.  The  concept  is 
nothing  more  than  the  specification  of  the  prefilter  and  the  feedback  compensa¬ 
tion,  depicted  in  Figure  1  of  Section  2.2,  in  terms  of  this  "optimum"  model. 

As  Campbell  points  out  (Reference  2),  these  "optimum"  characteristics  do  not 
peak  sharply,  and  reasonable  departures  are  allowable.  Thus,  the  approach 
seems  feasible. 

The  model  concept  as  applied  to  a  prefilter  is  illustrated  in  Figure  3. 


M 


-\ 

KG 

J — 

FIGURE  3  THE  MODEL  IN  THE  PREFILTER 


The  response  is  given  as 


and  if  is  sufficiently  high 


(2.12) 


(2.13) 
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The  model  concept  as  applied  to  feedback  compensation  is  illustrated 
in  Figure  4. 


c 


1 

1 

M 

FIGURE  4  INVERSE-MODEL  IN  THE  FEEDBACK 
The  response  is  given  as 

c  1.*#/ 
and  if  K  is  sufficiently  high  i 


(2.14) 


*  +  ~FT 


R  ~ 


=  M 


(2.15) 


Combinations  of  prefilter  and  feedback  compensation  can  be  used  to  in¬ 
dependently  control  the  response  to  commands  and  the  response  to  disturbances  - 
say  gusts.  If  we  designate  Mp  as  the  prefilter  model,  and  M ^  as  the  feedback 
model,  then  the  scheme  can  be  depicted  as  follows  if  the  gust  disturbance  is 
considered  to  affect  the  response  directly.  4 


FIGURE  5  COMBINED  PREFILTER  AND  FEEDBACK  MODELS 
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The  independent  responses  are  given  as 


and 


If  K  is  sufficiently  high 


mpmh 


KG 

_R_  -  Hjl-  M* 

6  "  KG  i  ,  KG 

nH 


(2.16) 


(2.17) 


(2.  18) 


(2.19) 


where  M ^ 


is  the  command  model,  and 


is  the  gust  (disturbance)  model. 


The  problems  associated  with  obtaining  a  desired  gust  response  are 
manifest  from  Equations  2.  18  and  2. 19.  The  feedback  model,  My  ,  and  com¬ 
pensation  in  the  forward  loop,  K  &  ,  provide  the  means  for  tailoring  the  gust 
response.  In  terms  of  frequency  spectra,  if  R/$  is  to  be  minimized,  then 
fly  must  be  minimized  over  the  expected  frequency  range  of  (j  .  A  high 

value  of  gain,  j(  ...will  tend  to  minimize  the  effect  of  gust  disturbances. 


Thus,  in  theory,  it  is  possible  to  obtain  any  desired  response  to  a  dis¬ 
turbance  input  (  R/j>  by  means  of  compensation  in  the  forward  and  feedback 
loop 8  provided  the  loop  gain  is  high.  Similarly,  in  theory,  it  is 

possible  to  obtain  any  desired  response  to  a  command  input  (  R/C  )  by  means 
of  compensation  in  the  prefilter  and  the  feedback  loop  ( Mp  My)  provided  the 
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loop  gain  is  high.  In  practice,  stability  problems  limit  the  value  of  gain 
that  can  be  used;  practical  problems  of  mechanization  limit  the  range  of  de¬ 
sired  responses  that  can  actually  be  realized. 

Other  reasons  for  combining  prefilter  and  inverse -feedback  model  tech¬ 
niques  arise  from  practical  considerations.  It  may  be  convenient  to  generate 
part  of  the  model  in  the  feedback  and  the  remainder  in  the  prefilter.  Or,  the 
prefilter  may  be  needed  to  cancel  undesired  feedback  characteristics  which  are 
generated  necessarily  for  purposes  other  than  inverse-model  generation.  Using 
both  types  of  compensation  will  generally  allow  more  flexibility. 

2. 4  .Aircraft  Model  and  Variables 

The  longitudinal  control  problem  was  picked  for  study  because 
difficulties  have  been  encountered  with  longitudinal  motions  in  the  past  (un¬ 
stable  airplane -pilot  combinations  have  been  encountered,  commonly  called 
pilot-induced  oscillations).  Also,  the  longitudinal  short-period  mode  is  one  of 
die  primary  control  modes  and  considerable  research  has  been  directed  toward 
determining  optimum  response  characteristics  for  this  mode.  A  considerable 
portion  of  this  research  has  been  performed  by  the  Cornell  Aeronautical  Labora¬ 
tory  and  is  reported  in  References  10  through  15.  The  earlier  results  are 
summarized  in  the  Appendix  of  Reference  2;  the  later  results  are  presented  in 
References  14  and  15.  The  results  of  most  of  the  research  on  longitudinal 
handling  qualities  done  to  date  (including  that  done  by  NASA,  various  govern¬ 
ment  agencies,  research  laboratories,  and  universities)  are  summarized  in  Ref¬ 
erence  16,  and  a  less  comprehensive  summary  is  found  in  Reference  17. 

The  data  of  Reference  14  define  combinations  of  short-period  undamped 
natural  frequency  and  damping  ratio  which  are  acceptable  to  the  pilot.  The  best 
tested  area  might  be  described  as  a  "plateau"  with  gently  sloping  sides  and  no 
distinct  "peaks"  or  optimum  combinations  of  CDn  and  £  .  In  selecting  a  model 
however,  consideration  must  be  given  to  the  relation  of  the  model's  speed  of 
response  or  bandwidth  with  respect  to  the  speed  of  response  or  bandwidth  of 
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the  closed-loop  system  incorporating  the  model.  The  bandwidth  of  the  closed- 
loop  system  must  be  at  least  as  large  as  that  of  the  model.  Particularly  if 
the  model  is  used  as  a  prefilter  and  it  is  desired  that  the  over-all  response 
be  that  of  the  model,  then  the  bandwidth  of  the  closed  loop  following  the  model 
(see  Figure  3  and  discussion  of  Section  2.3)  must  be  several  times  that  of  the 
model  alone.  Thus,  it  is  obvious  that  the  lower  the  bandwidth  of  the  model, 
the  less  difficult  will  be  the  design  of  the  system. 

The  short-period  handling  qualities  data  of  References  14  and  IS  are 
presented  in  Figure  10  in  the  form  of  a  complex-plane  plot.  This  figure  pre¬ 
sents  the  data  for  short*period  configurations  with  complex  roots  in  a  straight¬ 
forward  and  useful  manner.  For  configurations  with  real  roots,  the  plot  de¬ 
generates  to  a  single  line,  the  real  axis,  and  for  this  reason  this  representa¬ 
tion  was  not  used  in  References  14  or  15.  However,  for  the  purpose  of  se¬ 
lecting  an  adaptive  system  model,  it  is  of  some  value,  especially  when  root 
locus  techniques  are  used  for  synthesis.  Figure  10  shows  that  the  point  u)n 
a  1.0  cps  and  £  =  .9  is  in  the  "Best  Tested"  area  and  is  acceptable  from  the 
pilot's  point  of  view  as  a  short-period  model.  However,  the  point  &on  =  0.5 
cps  and  {  .7  is  also  in  or  on  the  boundary  of  the  "Best  Tested"  area  and, 

from  the  point  oft  view  of  the  adaptive  control  system  designer,  would  be  pre¬ 
ferred  because  of  its  lower  bandwidth. 


The  model  or  desired  response  that  results  from  this  research  may 
be  stated  in  transfer  function  form  as: 


n 

r 


An 


l  , 


K, 


~  incremental  normal  acceleration,  g 
=  stick  force,  lb 

=  short-period  frequency,  rad/sec 
=  snort-period  damping  ratio 
=  ( ^ | »  the  inverse  of  the  steady-state 
stick  f^rce  per  g. 


(2.20) 
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The  value  of  K-  depends  on  the  configuration  of  both  the  aircraft  and 
the  pilot's  control  stick.  For  conventional  lighter  airplanes,  with  high  maxi¬ 
mum  load  factors  (7  g  maximum),  typically  good  values  oi(Fs/n lie  be¬ 
tween  4  and  8  lb/g.  For  conventional  bombers,  with  low  permissible  load  fac¬ 
tors  (4  g  maximum),  typically  good  values  of  Jri/n^  lie  between  40  and  80  lb/g. 
It  appears  that  a  suitable  value  ot(Fs/t\^)Si  is  obtained  by  matching  the  max¬ 
imum  normal  acceleration  capabilities  of  the  airplane  to  the  maximum  control 
effort  expected  of  the  pilot. 


The  important  aspects  of  the  model  defined  by  Equation  2.20  are  listed 
for  emphasis. 


The  model  is  a  normal  acceleration  (  r>j)  one. 

The  frequency  (u>n  )  and  damping  [  'C  )  are  specified  constants 
invariant  of  the  airplane  configuration  or  flight  condition. 
However,  a  winged  vehicle  is  implied  which  derives  its  maneu¬ 
vering  capabilities  from  aerodynamic  lift  which  is  in  turn 
controlled  by  the  vehicle's  angle  of  attack. 

The  numerator  of  the  ^/Fs  transfer  function  -  more  generally 
n  /pilot  input  -  involves  no. time  constant  or  zeros. 

The  gain  of  the  A^/F%  transfer  function  is  constant,  though  it 
may  vary  from  one  vehicle  type  to  the  next. 


The  results  of  the  above  described  handling  qualities  research  have 
received  wide  distribution.  People  working  in  the  control  system  field,  par¬ 
ticularly  adaptive  control,  have  applied  these  results  to  the  model-reference 
concept.  The  particular  characteristics  for  the  transfer  function  have 

been  applied  to  all  manners  of  control  system  -  pitch  rate  and  roll  rate  com¬ 
mand  systems,  pitch  attitude  and  roll  attitude  command  systems,  etc  -  see 
Reference  3  for  examples.  This  ^.discriminate  use  of  the  n^/Fs  model  - 
while  not  objectionable  when  used  purely  for  illustrative  purposes  —  becomes 
objectionable  and  damaging  when  it  occurs  in  actual  flight  control  systems. 
Misuse  of  the  n^/F^  model  has  occurred  in  several  cases  where  it  has  been 
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applied  to  the  design  of  adaptive  pitch  rate  command  systems  -  &/ F^ 

systems.  When  evaluation  pilots  have  flown  these  systems,  they  have  com¬ 
plained  of  sluggish  response,  and  that  the  response  was  not  invariant  with 
flight  condition.  One  answer  to  complaints  of  sluggish  response  has  been  to 
reduce  the  damping  ratio,  to  as  low  as  f  =  .2  or  .3,  leading  to  the  conclu¬ 
sion  that  the  desirable  values  determined  by  the  handling  qualities  research 
are  in  error.  The  following  development  gives  a  better  appreciation  of  what 
is  involved,  and  why  care  must  be  used  in  selecting  the  reference  model. 

Three  responses  of  the  X-15  airplane  to  step  inputs  are  compared  in 
Figure  11:  one  is  for  the  uncontrolled  airplane;  the  other  two  are  for  the 
airplane  with  a  high -gain  control  system  which  employs  a  model  to  tailor  the 

response.  The  model  in  both  cases  is  the  )Tu  model  of  Equation  2.20  with 

rf  • 

=  .7;  but  in  one  case  it  is  used  as  a  Q  model,  in  the 

other  case  as  an  o£  model  (equivalent  to  an  model  if  Lg  is  small  and  the 

pilot  is  near  the  center  of  gravity).  It  is  clear  that  in  both  controlled  cases 

the  design  response  has  been  obtained;  in  one  case  the  model  response  is  dup- 

licated  by  ©  ,  in  the  other  by  Qt .  However,  to  one  familiar  with  handling 

qualities  requirements,  the  response  with  the  0  model  is  clearly  objectionable 

while  that  with  the  Q/C  model  is  highly  desirable.  The  objectionable  quality  in 

the  one  response  is  the  very  sluggish  response  in  ,  characterized  by  the 

very  long  time  constant  apparent  in  oC  . 

The  situation,  as  illustrated  in  Figure  11,  can  be  clarified  by  examining 
the  pertinent  block  diagrams  and  equations.  First  consider  the  "  Q  model"  case. 


e. 


Fs  50  9c 


AIRCRAFT  AND 
CONTROLLER 


■voC 


e 


FIGURE  6  BLOCK  DIAGRAM  WITH  INVERSE  0  -MODEL 

IN  T  IE  FEEDBACK 
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&  Response  (Figure  6): 

For  high  loop  gain 


&(s)  ~  jvj 

4.(0 ' 


e 


+  I 


where  U>n  -  .  5  cps 

r  =-7 


Response  (Figure  6): 

0C(i) 

A  (t) 

c 


cCCi) 

sJD  x  &u) 

6(*)  o  ($) 

C 


(2.21) 


(2.22) 


Substituting  the  airplane  transfer  functions,  given  by  Equations  A- 2  and  A- 3 
of  Appendix  A,  into  Equation  (2.  22),  yields 


j  +  0  &(s)  (2.23) 

~3C  (*)  Kfr  (T^  s+l)  6c  (s) 

Substituting  the  expressions  for  K ^  and  ,  given  by  Equations  A-l4  and  A- 13 
of  Appendix  A,  into  Equation  2.23,  yields 


<£M  /MS  *  MA  \ 


0(s) 


(2.24) 


where  j  and  '*'ot  are  given  by  the  following  expressions  (Equations  A- 14  and 
A- 15  from  Appendix  A). 


-L, 


Mr 


(2.25) 


*SL*-I\*LS 


(2.26) 


21 


ASD-TR-6 1-104,  Vol.  I 


The  relations  may  be  cimplified  by  neglecting  the  often  small  effect  of 
elevator  lift,  Lg  ,  to  obtain  the  following  approximation.  • 


*(  1  \  L  eC$) 
4  &  (Lbc)(T£+l)  $6.(s) 


Now  introducing  the  approximation  given  by  Equation  2.21,  yields 


CCU) 

ecis) 


(2.28) 


It  is  apparent  from  Equation  2.28  that  the  C/  response  will  be  dominated  by 
a  long  time  constant  approximately  equal  to  .  Also,  since  Lt*  changes 

with  flight  condition,  it  is  apparent  that  the  02  or  response  wiil  no.:  be 

invariant  with  flight  condition.  Furthermore,  the  M  cC  -model”  responses  of 

0 

Figure  1 1  show  that  a  pilot  does  desire  a  large  initial  overshoot  ir.  Q  ,  as  this 

is  necessary  in  order  to  give  the  desired  oC  and  6  response-.  If  the  pilot 

were  faced  with  the  response  given  using  the  -model  (Equation  2.20)  as  a 

Q  -model,  then  he  would  perhaps  settle  for  a  lower  damping  ratio  than  f  ST  .  7 

• 

in  order  to  obtain  this  initial  &  overshoot.  The  complaints  of  the  evaluation 
pilots,  referred  to  earlier,  are  thus  explained  and  can  be  attributed  to  improper 
application  of  the  data  obtained  from  the  handling  qualities  research. 

Next  consider  the  ”  -model”  case.  The  model,  in  this  case,  is  formed 
by  adding  a  prefilter  to  the  arrangement  of  Figure  6.  The  purpose  of  the  pre¬ 
filter  is  to  add  a  zero  or  a  numerator  time  constant  in  Qti)/&  £*)  so  as  to  cause 
the  proper  initial  overshoot  in  the  Q  response  to  a  step  command  input. 

* 
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FIGURE  7  BLOCK  DIAGRAM  WITH  INVERSE  Q  -DENOMINATOR 
IN  THE  FEEDBACK  AND  &  -NUMERATOR  IN  THE  PREFILTER 

# 

Q  Response  (Figure  7): 

a(s)  a(s)  6p(s)  (r;  $+■/)  ecs) 

4  (0  ^  6p(>)  ec(s)% (Tz  S+  l)  SpCs)  (2.29) 

For  high  loop  gain,  &/ & p  is  given  by  Equation  2.  21,  and  if  the  assigned 
values  for  2^  and  T"z  are  introduced,  then  Equation  2.  29  becomes 

Ahl-s  (n‘+'>  m- 

4  m  6  <2-30> 


Again  simplifying  by  neglecting  Lf  and  the  small  time  constant  Tz, 
the  following  expression  is  obtained 
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CL  Response  (Figure  71: 


The  &£  response  can  be  written  from  Equations  2.  23  and  2.29  as 
cC(  i)  (T,  $  +  !)(  s  l)  S(i) 

4  <■*>  *  K £  Crzs  -hOCr#  s  +  i)  6pU)  (2.32) 

If  the  assigned  values  for  t,  and  tz  are  introduced,  then  Equation  2.32 
becomes 


QC(i)  ££s) 

4(0  "  Kk  Spb)  1 

For  high  gain  6/&p  is  given  by  Equation  2.21,  and  if  Lg  is  again 
neglected.  Equation  2.  33  becomes 


\  l~cC  / 


(2.34) 


It  is  apparent  from  Equations  2.31  and  2.34  that  the  desired  results  are  achieved 
by  the  arrangement  depicted  in  Figure  7.  The  proper  Q  -model  numerator 
(Equation  2.31)  is  provided  to  give  the  proper  OC  -model  (Equation  2.34)  and 
hence  the  desired  -model  (Equation  2.2C).  Proper  control  gearing  is  pre¬ 
sumed,  and  the  desired  response  will  be  realized  only  if  the  loop  g  ’in  is 
sufficiently  high. 


In  summary,  both  systems  depicted  in  Figures  6  and  7  are'  pitch  rate 
command  systems,  but  essentially  different  reference  models  are  used.  The 
two  models  may  be  expressed,  without  gain  specification,  as: 


I 


ZT 


+ ) 


X  *+! 

W^TTT 


UNDESIRABLE 

DESIRABLE 


(2.35) 

(2.36) 
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The  desirability  of  the  two  models  is  evaluated  on  the  basis  that  the  input  comes 
from  the  pilot,  as  from  a  control  stick,  and  that  the  pilot  is  controlling  the  air¬ 
craft's  primary  motions  -  specifically  the  longitudinal  short-period  motions. 
Another  way  to  state  this  is  that  the  pilot  is  directly  in  the  micro-control  loop. 
For  other  inputs,  say  from  an  automatic  macro-controller  or  autopilot,  a  dif¬ 
ferent  model  might  be  desirable. 


The  above  results  may  be  related  more  closely  to  the  data  obtained  from 
the  handling  qualities  research,  and  the  story  made  complete  by  considering  the 
steady-state  (gain)  relations. 


A.  -  Ht  iZh.  *  A 

Fj  t  0/Se.  \  Fs 


(2, 37) 


Introducing  the  steady-state  values  for  &/Sg  ,  ir/St-  .  and  y*  from  Equations 

A-7,  A-9»  and  A-lOof  Appendix  A,  and  for  *l/ps  from  Equation  2.20,  and  letting 
OjOc  =  1 ,  then  the  recommended  gearing  (gain)  is  found  to  be 


-  _2_  k 

V  Knj- 


(2.38) 


where  K_  ,  the  inverse  of  the  steady-state  stick  force  per  g,  is  an  appropri- 
ately  selected  constant  as  discussed  earlier. 


The  foregoing  analysis  has  all  been  based  on  the  assumption  that  die  de- 
sired  response  for  Q  can  be  achieved  (i.  e. ,  that  the  aircraft's  response  will 
duplicate  that  of  the  reference  model).  Even  with  this  assumption,  it  is  neces¬ 
sary  to  program  the  control  gearing  (  )  with  speed  (  V ),  and  the  numerator 

time  constant  of  the  model  (  7^  )  with  flight  condition  on  the  following  basis. 


(2.39) 


25 


ASD-TR-6 1-104,  Vol.  I 


The  gearing  can  be  programmed  with  speed  rather  easily,  and  accuracy  is  not 

important.  But  programming  the  time  constant  in  the  manner  indicated  is 

difficult  and  objectionable  (particularly  if  the  approximation  of  7$  by  (A*  cannot 

be  allowed).  Furthermore,  in  practical  applications  the  ideal  response 

0 

(model  response)  for  &  will  not  usually  be  achieved  due  to  gain  limitations, 
and  the  required  values  for  the  time  constant  (  )  may  thereby  be  affected. 

(For  example,  a  specific  case  rvbere  this  occurs  is  dealt  with  in  Section  3.3, 
where  the  G. E.  System  is  applied  to  the  X-15.)  Accordingly,  a  better  method 
for  obtaining  the  desired  C>C  or  fiy  response  may  be  to  use  dC  or  feed¬ 
back  directly  in  the  control  system.  This  is  an  appropriate  subject  for  further 
study. 


The  control  problem,  for  the  longitudinal  short-period  motions,  can  be 
formulated  in  rather  genera)  fashion.  One  can  categorize  the  types  of  command 
inputs,  the  variables  that  are  suitable  for  sensing,  the  variables  that  need  to  be 
controlled  in  some  satisfactory  or  desired  manner,  the  suitable  control  quan¬ 
tities,  and  the  tasks  or  objectives. 


I 

Types  of  Command  Inputs 

a. 

Pilot  control 

1.  Force 

2.  Position 

b. 

Automatic  Control 

c. 

Monitored  Control 

11 

Suitable  Responses  for  Sensing 

a. 

Angle  of  attack, 

b. 

Normal  acceleration, 

c. 

Pitch  rate, 

d. 

Derivatives  of  (a),  (b),  and  (c) 

e. 

Integrals  of  (a),  (b),  and  (c) 

III 

Controlled  Responses 

a. 

Angle  of  attack,  o£ 

b. 

Normal  acceleration, 
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c.  Pitch  rate,  ^ 

d.  Flight  path  curvature,  if 

e.  Attitude  angle,  Q 

f.  Flight  path  angle,  Tf 

g.  Altitude,  ti 

IV  Suitable  Control  Quantities 

a.  Elevator  angle,  Sg 

b.  Reaction  control 

V  Task  or  Objective  -  Vehicle  Control  During: 

a.  Boost 

b.  Glide  or  Enroute  Flight 

c.  Re-entry 

d.  Approach 

e.  Landing 


Only  one  control  input  is  available  for  the  control  problem  under  consideration, 
so  that  only  one  possible  form  of  response  is  available  during  any  particular 
time  interval  for  all  of  the  variables  listed  in  Item  in.  Thus,  the  control  synthesis 
problem  can  be  reduced  to  the  following  steps:  select  a  suitable  model  for  the 
response  of  any  one  variable  (listed  in  in)  to  the  given  type  of  command  input 
(listed  in  I);  sense  one  or  more  of  the  variables  listed  in  II;  finally,  develop 
a  control  system  which  links  the  sensor  ( II  )  and  the  input  to  the  airframe  con¬ 
troller  (elevator  and/or  reaction  control),  and  use  the  model  as  the  reference 
for  the  desired  response.  The  task  or  objectives  are  influential  in  the  process 
of  selecting  die  model,  the  specific  application  governs  what  sensors  are  used, 
while  the  controller  configuration  is  dictated  primarily  by  the  performance 
required  of  the  system. 
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SECTION  III 

INITIAL  STUDY  AND  APPLICATION  OF  G.E.  SYSTEM 

3.  1  Introduction 

The  analysis  of  the  G.E.  system  was  conducted  in  two  separate  phases. 

The  phase  reported  in  this  section  ( Section  III)  was  a  general  study  of  the  sys¬ 
tem  and  its  application.  This  study  used  linear  analysis  techniques  based  prin¬ 
cipally  on  the  root  locus  method.  An  ESLAC  computer  was  used  to  obtain  the 
root  locus  plots.  This  computer  and  the  form  of  the  root  locus  plots  are  des¬ 
cribed  in  Appendix  B.  The  system  also  was  simulated  on  an  analog  computer. 

The  analog  computer  used  to  simulate  the  system  in  the  first  study  was  a  small 
EASE  computer.  Simulation  of  the  airplane,  the  actuator,  the  reference  model, 
and  the  frequency  sensor  used  its  full  capacity  and  few  components  were  avail¬ 
able  for  simulating  any  additional  elements  or  nonlinearities. 

The  second  phase  of  the  analysis  was  directed  at  specific  problems  in  ap¬ 
plication  of  the  G.E.  system.  The  results  of  this  study  are  reported  in  Section 
IV.  The  ESLAC  computer  was  again  used  in  the  analysis  and  the  system  was 
simulated  on  a  larger  REAC  analog  computer  which  permitted  simulation  of 
actuator  nonlinearitiea. 

The  X-15  was  the  study  vehicle  in  both  phases.  The  G.E.  system  was  ap¬ 
plied  to  the  problem  of  controlling  iae  X-15  short-period  dynamics  during  re-entry. 
The  equations  used  to  represent  the  airplane  are  contained  in  Appendix  A. 

During  the  first  study  using  the  ESLAC  and  EASE  analog  computer  the 
following  problem  areas  were  examined. 

1.  General  Analysis  of  the  System 

a.  Block  diagrams 

b.  Model  concept 

c.  Root  locus  analysis 

2.  Application  of  the  G.E.  system  to  the  X-15 

a.  Selection  of' system  parameters 

b.  Response  to  step  command  inputs 

c.  Response  to  gusts 

d.  Modification  of  basic  airplane  characteristics  -  particularly 
unstable  configurations 
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e.  Sensor  and  actuator  dynamics 

3.  Adaptive  loop  characteristics. 

During  the  second  study  using  the  ESXAC  and  the  REAC  analog  com¬ 
puter  the  following  specific  problems  were  examined. 

1.  The  response  of  the  system  to  inputs  at  the  actuator  valve 

2.  The  effect  of  actuator  nonlinearities 

3.  Characteristics  of  the  frequency  sensor 

4.  The  effect  of  "noise"  in  the  rate  gyro  output 

3.2  General  Analysis  of  the  System 

The  General  Electric  Self-Adaptive  Flight  Control  System  (GESAC)  - 
referred  to  in  this  study  as  die  G.E.  System  -  involves  the  following  general 
concepts: 

1.  Inverse  model  in  the  feedback, 

2.  Rate  gyro  sensing  device, 

3.  Invariant  &  response  with  flight  condition, 

4.  Adaptive  control  of  the  primary  or  forward  loop 
gain  employing  a  frequency  sensor,  and 

5.  No  inserted  test  signal,  per  se. 

Two  block  diagrams  of  the  system,  applied  as  a  pitch  rate  control  system,  are 
presented  in  Figure  8.  One  shows  the  system  as  actually  mechanized,  the  other 
•hows  die  system  as  it  would  appear  with  the  model  transfer  function  explicitly 
in  the  feedback.  The  purpose  of  the  K^/s  block  in  the  forward  loop  is  to  pro¬ 
vide  a  zero  error  system,  so  the  gain  of  the  closed-loop  system  will  be  invariant 
of  the  aircraft  gain. 


The  equations  for  the  model,  in  terms  of  the  parameters  of  the  physical 
•ystem  depicted  in  Figure  8,  are 


i* 


' K+Tc  +KZTC  ' 

'  K,TC  + 

*Sl  K,  J 

*i 

m 

(3.1) 

*TC 

+  / 

Uz)  =  K‘ 

(3.2) 

M  T&  (K++Kz) 


29 


ID-1471-F-1 


ASD-TR-b  1-104,  Vol.  I 


K/  Tc  +  Kz 

K, 


(3.3) 


There  are  four  parameters  available  ( 
theory,  K,  ,TC  •  and 

there  are  practical  limits. 


K.l  ,  Kz,  ,  and  Tf,  ),  so  that  in 
can  all  be  set  independently  -  actually 


The  general  concept  of  how  the  system  works  is  based  on  the  explana¬ 
tion  given  in  Section  2.  3  j  the  gain  K3  is  kept  high  enough  so  that  the  system 
response  approximates  that  of  the  model;  the  function  of  the  frequency 
sensor  is  to  maintain  the  gain  at  an  optimum  value  -  as  high  as  possible  without 
causing  instability.  The  specific  operation  of  the  system  can  best  be  described 
using  the  root  locus  technique.  Two  root  locus  plots  are  presented  in  Figure  12, 
one  for  a  second-order  actuator,  the  other  for  a  first-order  actuator.  To  illus¬ 
trate  the  general  concept  only  the  actuator,  airplane,  and  model  zeros  and  poles 
need  be  included. 

Consider  first  the  case  with  the  second-order  actuator  as  depicted  in 
Figure  12.  As  the  loop  gain  of  the  system  is  increased,  the  aircraft  poles  ap¬ 
proach  the  feedback  zeros  (the  poles  of  the  model),  the  integrator  pole  at  the 
origin  approaches  the  airplane  zero  (  at  ),  the  feedback  pole  (at  -//7j  ) 

goes  out  the  negative  real  axis,  and  the  actuator  poles  approach  the  imaginary 
axis  and  eventually  go  into  the  right-half  plane  -  unstable.  The  limiting  value 
of  the  loop  gain.is  determined  by  the  actuator  poles,  and  some  value  less  than 
that  for  instability  represents  the  optimum  value  of  the  gain. 

In  order  to  see  how  the  system  works  and  understand  the  design  phil¬ 
osophy,  it  is  necessary  to  consider  the  equations  which  define  the  open-loop 
and  closed-loop  response  in  terms  of  the  zeros  and  poles  of  the  root  locus  plot 
(Figure  12.,  also  refer  to  Figure  8). 
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Forward  Loop: 


■£%  .L«fa4f 

*6.)  F  *  ¥  Dp(*> 


(3.4) 


Kf  -  Kt  K3  Kmt  Ka/c 


(3.5) 


Npk)  =  I  +  ,/yJ  j 


°f(0-  5 


»/*S 


Feedback  Loop: 


(4?+^Fs+j) 

4crv  ja/c 


(3.6) 


(3.7) 


ac'sj/  (s)sMl 

6(s)  KHbH^ 


(3.8) 


K„=l 


(3.9) 


co, 


Is.,) 

»  /t 


MODCL 


(3.  10) 


Cloaed  Loop: 


Dhw  =  (l  +  ifa  ) 


^6)  Kp&p(s)  Kf  Mf  (s)  Dh  ($) 


(3.11) 


1+ KfGf (sJKhChC*)"  0HCs)DF(s)^KfNp(ONH(s) 


(3.12) 


(3.13) 
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[ 

1 

The  root  locus  plot  determines  the  zeros  of  OCi)  which  are  the  poles  of  (?(}).  I 

The  zeros  of  £(s)  are  given  directly  by  the  zeros  of  Np( s)  and  .  The 

closed-loop  gain  is  obtained  easily  from  the  zero  frequency  value  of  the  closed- 
loop  response: 


Li  m 

S-*0  Oc(s)  (3.14) 

Thus  the  following  closed-loop  characteristics  are  obtained: 


Lim 

S-+0 


by  Cs)  bpCs)  4.  Kf  Wp  (i)  Ntf  (s) 


5  —>0  $ +Kp 


K  *  / 


(3.15) 

(3.16) 


NC0-(l*i fo)(l+lfc) 


(3.17) 


=  pole  in  region  O 
=  magnitude  of  the  pole  in  Area  A 
=  measure  of  the  angle  of  the  pole  in  Area  A 
=  magnitude  of  pole  in  Area  B 
=  measure  of  the  angle  of  the  pole  in  Area  B 
=  pole  along  the  real  axis  emanating  from 


(3.  18) 


The  design  approach  assumes  that  the  poles  emanating  from  the  actuator 
can  be  maintained  in  Area  B  by  appropriately  adjusting  the  loop  gain,  and  that 
correspondingly  the  following  conditions  will  be  obtained: 
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i. 

(3.19) 

L 

'  *  c 

(3.20) 

»  /jj 

B  "ft 

(3.21) 
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so  that 


a  ~ 


Tc  S+I 


* 


TT 


s  -hi 


(3.22) 


The  desired  response  is  that  given  by  Equation  3.22,  and  the  question  is,  how 
stringently  must  Equations  3. 19  -  3.21  hold  in  order  to  achieve  this?  Truxal 
(Reference  9,  p.  43)  gives  the  criteria  that  if  is  to  have  negligible  effect 
on  the  response,  then 


Using  the  same  reasoning,  we  may  say  that 

"»e  >  14 


(3.23) 


(3.24) 


The  criteria  on  ^  poses  a  much  more  significant  problem.  Since  this  pole 
is  near  die  origin,  it  could  dominate  the  response  if  it  were  not  close  enough  to 
the  zero  at  ,  and  thereby  make  the  adaptive  system  ineffective.  It 

is  necessary  to  evaluate  the  residue  at  this  pole  in  order  to  determine  its  effect. 
The  residue  for  a  unit  step  response  is 


R 


=■  [(s+?>0)i-  GU)]S  s  _ 


ro 


(3.25) 
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The  factors  in  G(s)  associated  with  the  poles  and  zeros  that  are  far  from  the 
origin  will  contribute  negligibly  to  R_ 


5 


>/-Tq  )(l+  t/rt,  ) 


/  s* 

zrt 

\ 

5 

5 

± .. ,  \f 

i  - 

To 

5= 


*0 


(3,  26) 


(3.27) 


To  obtain  some  feel  for  the  numerical  values,  consider  a  case  where  ,  y  /  , 

and  are  numerically  small  so  thac 


n  «  A  , 

S  -  775  ■' 


-  /  - 


.  JL- 


'/*! 


(3.28) 


Since  that  part  of  the  amplitude  of  the  step  response  associated  with  the  pole 
is  just  f?  -  ,  then  if  -70  0  <1/2  ( tke  8teP  response  will  be  dominated 

'O  j  f 

by  the  pole  at  -jo  ,  If  ^  ,  then  perhaps  should  not  differ 

**  ^c, 

from  ///%  by  more  than  10  or  20%.  As  decreases,  then  this  percentage 
difference  should  be  even  less.  Thus,  it  is  very  important  that  the  pole 
at  should  be  near  the  zero  at  .  Hence  Equation  3. 19  is  a  strong 

requirement. 


The  response  of  the  closed-loop  system  is  characterized  by: 

1.  Short-period  mode  defined  by  the  roots  found  in  area  C. 

2.  Actuator  mode  defined  by  the  roots  found  in  area.  B. 

3.  A  numerator  time  constant,  T0  --  which  actually  turns  out 
to  be  quite  small  and  of  negligible  effect. 

4.  A  pole  and  a  zero  in  region  D. 


34 


ASD-TR-6 1-104,  Vol.  I 


'S 


5.  Unity  closed-loop  gain.  It  the  feedback  gain,  ,  is  not  unity, 
then  the  closed-loop  gain  would  be  >/ . 

The  function  of  the  adaptive  device,  the  frequency  sensor,  is  to  main¬ 
tain  the  gain  at  the  optimum  value  -  in  area  B.  The  frequency  sensor  measures 
the  response  of  the  closed-loop  actuator  mode,  compares  it  with  a  reference 
frequency  O>0  ,  and  changes  the  loop  gain  so  as  to  null  the  difference  in  the  two 
frequencies.  This  can  be  expressed  symbolically  as  follows: 

s  4  (to I  -6J.) 

dt  ‘  1  °J  (3.29) 


where  cut  is  the  frequency  of  the  input  signal  to  the  frequency  sensor.  If  the 
input  frequency  is  too  high,  the  gain  decreases;  if  the  input  frequency  is  too 
low,  the  gain  increases.  The  value  of  the  fiequency  is  picked  so  that  the 
damping  of  the  actuator  mode  will  be  within  the  limits 

.2<  r<.3  (3.30) 

The  actuator  mode  is  not  supposed  to  go  unstable,  but  on  the  other  hand,  u  def¬ 
inite  oscillation  is  required  in  order  to  sense  the  frequency.  The  details  of 
the  frequency  sensor  mechanization  are  d  cussed  in  Sections  3.4  and  4.4. 

The  significance  of  the  four  areas  in  Figure  12  is  as  follows: 

Ar«*a  ."A":  This  is  the  envelope  of  the  aircraft's  short*period 

roots  throughout' the  flight  spectrum. 

Area  "B":  This  is  the  envelope  in  which  the  closed-loop  actuator 

mode  is  maintained  by  the  frequency  sensor.  Changes 
in  frequency  are  associated  with  noine  and  frequency 
sensor  resolution.  Changes  in  damping  are  associated 
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with  changes  in  the  actuator  locus  caused  by 
changes  in  the  aircraft  poles  or  other  poles 
and  zeros. 

Area  "C":  This  is  the  envelope  for  the  short*period  roots  of 

this  controlled  airplane.  This  envelope,  then,  rep¬ 
resents  the  actual  adaptive  capabilities  of  the 
system.  The  feedback  zeros  may  be  thought  of  as 
the  model  generators  rather  than  the  true  model. 

The  system  is  designed  so  that  the  desired  model 
roots  fall  in  area  C,  and  if  area  C  can  be  made  small 
then  the  system  provides  invariant  response. 

Region  "D":  This  is  the  region  or  range  for  the  pole-zero  com¬ 

bination.  This  pole -zero  combination  can  have  a 
significant  effect  on  the  closed-loop  transient  re¬ 
sponse  if  they  do  not  cancel.  How  well  they  cancel 
depends  on  the  maximum  gain  achieved.  This  pole- 
zero  combination  is  of  particular  significance  if  a 
prefilter  is  used  to  shape  the  response  as  described 
in  Section  2.4. 

The  second  c£te  presented  in  Figure  12  employs  a  first  order  power 
actuator  and  a  second  order  servo  to  control  the  power  actuator  -  similar  to 
the  X-15  control  system.  Modes  of  higher  frequency  than  the  actuator  need  to 
be  included  in  this  **ase  in  order  to  cause  the  "actuator  modes"  to  go  unstable  - 
a  representation  of  rate  gyro  dynamics  would  have  served  equally  well. 

The  two  cases  illustrated  in  Figure  12  have  substantially  different  types 
of  root  loci,  and  the  comparison  points  out  some  of  the  design  problems. 

Cleanly,  the  general  form  of  the  root  loci  must  be  known,  and  the  locus  of  the 
"actuator  mode"  in  the  vicinity  of  the  operation  point  must  be  carefully  estab¬ 
lished  and  stabilized  (not  allowed  to  move  much).  The  operation  of  the  frequency 


i 
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sensor  depends  on  having  a  sufficient  gradient  of  frequency  with  loop  gain 
in  the  vicinity  of  the  operating  point.  That  is,  must  be  large,  where 

U)t  is  the  closed-loop  actuator  mode  and  if  is  the  open-loop  gain.  The  root 
locus  approach,  as  described  here,  sheds  little  light  on  the  dynamic  charac¬ 
teristics  associated  with  the  adaptive  loop  itself,  the  loop  through  the  frequency 
sensor.  This  subject  is  considered  in  Sections  3.4  and  Sections  4.4  and  4.5. 

3.3  Application  of  the  G.  E.  System  to  the  X-15  Airplane 
3.  3. 1  Selection  of  Parameters 

The  first  problem  considered,  in  applying  the  G.E.  System  to  con¬ 
trolling  the  longitudinal  shorfc>period  motions  of  the  X-15  airplane,  was  the  syn¬ 
thesis  of  a  control  system  of  the  form  shown  in  the  block  diagram  of  Figure  8. 
The  following  transfer  functions  were  assumed  in  the  various  blocks. 


Aircraft: 


fa  s+i) 


s+i 


Parameters  listed  in  Table  in. 


Actuator: 


Rate  Gyro: 


s~ 


e 


-  / 


(A-7) 


(3.31) 


37 


ASD-TR-6 1-104,  Vol.  I 


The  feedback  gains  and  time  constants  were  selected  to  give  the  following 
equivalent  model  in  the  feedback. 


Model: 


f±\ _ •*+' 

'  ec)M  &?) *  -gf-  s  ~i 


(3.32) 


(^n)^  -  .5  cps  =  3.14  rad/Gec,  =  .  V 


Feedback  Gafins  and  Time  Constants: 


K  -i 

*z  *  ^  [(££)"  "Te  1 


& 


CO 

n 


=  Variable  Gain 


from  Equation  3.3  (3.33) 

from  Equation  3.2  (3.34) 


To  establish  some  idea  of  the  possible  values  of  Tc  ,  the  parameters  for 
the  X-15  at  t  =  90  sec  in  the  re-entry  (see  Tables  I,  II,  and  III)  were  selected, 
and  a  root  locus  plot  was  made  for  the  closed-loop  system  with  various  values 
of  T0  ,  shown  in  Figure  13.  The  coordinates  of  Figure  13  are  those  used  by 
the  ESIAC  computer,  and  are  described  in  Appendix  B.  The  results,  when  in¬ 
terpreted  with  the  adaptive  concept  in  mind,  show  that  the  possible  values  of  Tc 
are  restricted.  Particularly,  note  the  case  for  7^.  =  .  5  and  the  locus  for  the 
act'jator  pole.  As  discussed  in  Section  3.2  ,  the  closed-loop  actuator  pole  must 
lie  in  the  region 


,2  <r<.3 

and  in  this  region,  the  actuator  frequency  must  increase  with  gain  ( 

This  is  not  the  case  for  =  .5.  In  fact,  there  are  two  possible  values  of  gain 
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(  K,  =  1.5  and  Kj  <  1)  for  the  proper  reference  frequency^  "S'  26  rad/sec. 
Furthermore,  the  aircraft  pole  is  far  from  the  model  zero  even  at  Kj  -  2.5. 
Clearly  Tc  must  not  be  too  large  -  the  pole  at  should  lie  to  the  right  (in 

EASIAC  coordinates)  of  the  actuator  poles,  or  >  (u>n )Acr  . 

Based  on  the  foregoing  results,  the  following  values  of  control  system 
constants  were  selected  for  use  with  the  second-order  actuator. 

Tc  =  .025  sec 

K,  =  l 

Kz  =  .413 
K*  =  3.39 

3. 3. 2  Response  to  Step  Commands 

The  performance  of  the  control  system,  from  the  linear  standpoint,  is 
shown  by  the  root  locus  plots  in  Figure  14  for  the  full  gamut  of  the  X-15's  flight 
range.  The  dynamic  pressure,  ,  ranges  from  20  psf  at  the  apex  of  the  alti¬ 
tude  run  to  a  1000  psf  maximum,  giving  a  200:1  gain  change  requirement.  The 
range  for  the  roots  of  the  short-period  mode  of  the  controlled  airplane  is  found 
to  be  .  35  cps  <.  «*>„<  .65  cps,  and  .4  <  X <.5. 

The  response  realized  with  the  control  system  configuration  of  Figure 
8  is  shown  in  Figure  11,  the"  Q  -model"  case.  The  various  model  concepts 
have  been  discussed  in  Section  2.4  and  need  no  elaboration.  An  attempt  was 
made  to  set  7^  *  Tj  in  order  to  produce  the  response  labeled  "of  -model", 
but  with  no  success.  (7^  must  be  used  to  shape  the  locus  of  the  high  frequency 
modes  as  discussed  previously).  Also  a  lead-lag  filter  ^ was  con¬ 
sidered  for  the  feedback,  but  the  necessary  separation  between  'TJ  and  ^ 
caused  an  intolerable  gain  loss  when  used  in  this  manner.  Accordingly,  a 
prefilter  was  inserted  (as  suggested  in  Section  2.4,  Figure  7  and  Equation  2.29) 
ahead  c t  the  input  summing  point. 

,$p/±L  *  (3.35) 

$c(s) 


39 


ASD-TR-6 1-104,  Vol.  I 


The  time  constant  was  first  selected  on  the  basis  of  Equation  2.30.  and 
7^  was  made  small.  The  response  for  the  X-15  at  t  =  90  sec  was  computed 
on  the  analog  using  the  following  prefilter  time  constants. 


"Tf  *  Tg  -  4.4  sec 
7*  *  *6/  Tj  -  .04  sec 


Surprisingly,  the  o£  response,  though  improved,  still  exhibited  the  long  time 
constant.  Analysis  of  the  situation  showed  that  it  was  the  pole-zero  combina¬ 
tion  (in  Region  D  of  the  root  locus  plot  of  Figure  12)  in  the  &/$c  transfer 
function  that  was  the  culprit  -  the  two  were  not  cancelling.  To  show  the  effect, 
the  complete  transfer  function  for  $  must  be  used  and  not  just  the  model 
equation.  Continuing  with  the  notation  of  Figure  7,  the  Q  response  can  be  ex¬ 
pressed  for  the  bandwidth  of  interest  using  Equations  3. 13and3.  Has 


Sis) 


( S  -/•  I ) 


( 


•2  !, 


w  *  ^ 


**•)('+ ft) 


(3.36) 


and  repeating  Eqviation  2.  32,  but  neglecting  the  small  time  constants  and 


OC(s)  _  (T,  s+/)  9  Ci) 

Ki%)  ^3 s+/^  4^ 


Now  when  Equations  3.  36  and  3.  37  are  combined,  we  see  that  the  pole  at  -  l/^Q 
(in  Equation  3.37)  is  cancelled,  leaving  the  pole  at  ~yo^  as  the  only  pole  in  the 
vicinity. 


oca) 

OcCs) 


(3.38) 
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Clearly,  the  desired  response  can  be  achieved  only  if 

77  “  (3.39) 

The  situation  is  illustrated  in  the  root  locus  plot  of  Figure  17,  part  a,  where 
the  poles  corresponding  to  the  appropriate  gain  {  K 2  )  for  % Q  =  .  2  is  indi¬ 
cated.  For  this  case,  the  X-15  at  t  =  90  sec, 

i°*  =  .1 

//^  =.23 

and  if  Yj  is  set  equal  to  ,  then  ^  <  ;£^//77jand  (based  on  the  approxi¬ 
mate  relation  of  Equation  3.28)  more  than  half  the  magnitude  of  the  o£  response 
will  be  comprised  of  the  sluggish  response  yielded  by  .  Thus,  in  the  prac¬ 
tical  situation  illustrated  in  Figure  17,  part  a,  and  also  for  the  06  -model  case 
of  Figure  11,  the  prefilter  time  constant  was  made  =  10  sec  rather  than 
V  *  ’S  =  <-5  sec  as  it  would  have  been  had  the  loop  gain  been  very  high. 

An  idea  of  the  error  involved  by  making  *  7^  can  be  obtained  by  looking  at 
the  response  in  Figure  15-  for  t  =  60  sec.  For  this  condition  ^  =  10  sec  was 
used  and  7^  =  11  sec.  It  is  seen  that  there  is  a  substantial  "tail"  in  oC  . 

Thus,  in  order  to  achieve  the  desired  response  for  oC  and  O  ,  it  is  necessary 
to  adapt  77  according  to  Equation  3.39.  Since ^  is  a  function  not  only  of 
7^  but  also  of  the  variable  gain  Kj  ,  this  becomes  an  even  more  difficult  task 
than  adapting  equal  to  7^  *  f/ . 

Responses  are  shown  in  Figure  15,  comparing  the  open-loop  airplane 
with  the  controlled  airplane,  for  all  the  selected  re-entry  times.  The  prefilter 
time  constants  were  selected  to  give  the  desired  response  for  the  X-15  at 
t  =  90  sec. 

*7^  =  10  sec 
=0.1  sec 
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The  prefilter  (  Tf  and  Tt  )  was  not  varied  for  two  reasons:  one,  to  evaluate 
what  the  effect  of  a  constant  rather  than  an  adapted  prefilter  would  be;  two, 
so  that  the  capabilities  of  the  system  to  maintain  invariant  G  response  could 
be  evaluated.  The  responses  shown  in  Figure  15  correspond  to  the  root  loci 
of  Figure  14.  The  variations  in  the  6  response  are  small;  a  long  a:  response 
time  appears  for  the  higher  altitudes,  but  the  prefilter  is  effective  to  a  certain 
extent  even  without  adapting  the  time  constant. 

3.  3.  3  Response  to  Gusts 

The  gust  response  of  the  airplane  with  the  G.  E.  System  is  compared  in 
Figure  16  with  that  of  the  uncontrolled  airplane,  and  with  that  for  a  control 
system  employing  just  high-gain  pitch-rate  feedback.  The  difference  between 
the  two  types  of  feedback  was  discussed  in  Sections  Z.  Z  and  Z.  3,  and  they  were 
shown  to  differ  principally  in  the  response  to  disturbances.  For  gust  inputs  of 
the  type  shown  in  Figure  16,  there  seems  to  be  little  difference  in  the  response 
for  the  two  types  of  control  systems. 


3.  3.  4  Modification  of  Basic  Airplane  Characteristics  - 
Particularly  Unstable  Configurations 


Configurations  suitable  for  hypersonic  flight  tend  to  have  unstable  pitching 
moment  curves  (positive  CM  ),  particularly  at  low  speeds  and  at  high  angles 
of  attack.  Accordingly,  the  X- 15  equations  were  modified  by  arbitrarily  changing 
Me  in  the  direction  of  static  instability,  and  both  root  locus  plots  and  analog 
responses  were  obtained  for  these  conditions.  To  complete  the  picture,  and 

were  also  varied  to  produce  dynamic  instability.  The  results  are  presented 
in  Figure  17  and  show  that  the  control  system  stabilized  the  response  in  all 


cases.  The  &  response  was  kept  invariant,  qualitatively,  but  the  initial  Q 
overshoot  increased  for  =  0  or  >  0.  The  root  locus  plots  show  that  the 

pole-zero  configuration  changes  substantially  for  >  0  and  therefore,  the 
prefilter  zero  at  ~1/T  -  -0.  1  does  not  cancel  the  real  pole  that  is  approaching 
,  but  now  from  the  high  frequency  side.  This  accounts  for  the  increase 
in  the  initial  Q  overshoot  and  the  attendant  overshoot  in  CC . 
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3.3.5  Sensor  and  Actuator  Dynamics 

The  results  discussed  in  the  preceding  paragraphs  of  Section  3.3  were  all 
obtained  with  the  seconctorder  actuator  (Equation  3.  3-1}  and  no  sensor  dynamics* 
The  effect  of  incorporating  an  actuator  of  the  X-15  type  and  rate  gyro  dynamics 
were  studied  on  the  ESIAC. 


Actuator: 


*  / 

5*JV' 


Xcr 


power 

actuator 


_ ^iCRVO _ 

actuator  servo 


) 


K*C  T 


*HKvo 


*  S^RVO 


K  1 

=  200  rad/sec  =  32  cps 
=  .4 


.  5  sec  >X^ct  >  *  1  **c 


(3.40) 


The  change  in  die  actuator  time  constant  is  produced  by  nonlinearity  in  the 

control  system:  'T‘act  -  -5  for  small  deflections,  CT  =  .1  for  large  deflections. 
This  nonlinearity  was  simulated  in  the  REAC  analog  studies  and  is  discussed 

in  Section  4.  3>.  # 

Qa  4  ^  US 

■UteQyxo:  =  (3.41) 


=  1 


y 0 


100  rad /sec  =  16  cps 
.  5 


The  effect  of  incorporating  the  above  actuator  and  sensor  characteristics 
into  the  system  is  shown  in  fee  root  locus  plots  of  Figure  18.  One  important 
observation  is  that  it  is  now  fee  poles  at  //V^cr  and  ;/7^  that  combine  and 
go  unstable,  and  it  is  fee  frequency  of  this  mode  that  must  now  be  detected  by 
th«  frequency  sensor,  and  used  to  adjust  the  system  gain.  The  two  root  loci 
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for  r>fiT  =  .5  and  .  1  sec  are  compared  in  Figure  18  and  show  that  a  problem 
is  presented.  The  variation  in  *r^6r  causes  the  root  locus  of  the  detected  mode 
to  vary,  and  no  satisfactory  reference  frequency  (u>0  )  can  be  picked.  If  the 
reference  frequency  u>0  is  picked  on  the  basis  of  the  YacT  =0.5  sec  locus 
(  eop  •  40  rad/sec),  then  the  frequency  sensor  will  attempt  to  provide  the 
following  gains. 

K3  =  20  for 

Kj  =  4  for 

The  aircraft  poles  are  close  to  the  model  zeros  for  K 3  -  20  but  not  so  close 
for  Kj  =  4.  The  given  situation  is  marginal,  but  if  the  variation  were 

larger  (say  .  5  <r^CT  <.05)  then  the  aircraft  response  would  be  hardly  adapted  at 
all  for  the  low  gain  case. 

An  attempt  was  made  to  compensate  for  the  Kcr  changes  by  inserting 
lead-lag  compensation  (of  the  form  )  in  the  control  system,  and  the 

results  are  presented  in  Figure  19.  r?he  lead  zero  was  centered  in  the  range  of 
the  actuator  pole  (  //y£er  )  so  as  to  attract  it.  The  lag  pole  was  placed  to  the 
left  (ESIAC  coordinates)  of  the  feedback  pole  (  >/Tc,  )  to  attract  it.  The  purpose 
was  to  stabilize  the  locus  of  the  mode  detected  by  the  frequency  sensor.  Fig¬ 
ure  19  shows  that  stabilization  was  accomplished  in  a  sense  -  the  two  loci  are 
much  closer  together  than  in  Figure  21.  Now,  however,  the  shape  of  the  loci  of 
the  detected  mode  is  poor  since  the  frequency  variation  with  gain  («£«£'  / at  K  ) 
is  small;  but  this  could  be  cured  by  introducing  an  extremely  fast  rate  gyro  - 
say  one  with  St  400  rad/ sec. 

However,  another  problem  exists  (demonstrated  in  both  Figures  18 
and  19)  which  may  impose  a  serious  limitation.  Consider  Figure  18  and  the 
conditions  imposed  (Equations  3.42)  by  selecting  -  40  rad/sec.  Assume 
that  die  system  is  operating  at  small  amplitudes  (no  pilot  input,  and  only  noise 
is  exciting  the  system)  so  that  it  would  be  operating  on  the  locus  for  7^^-  .  5  sec 


"'/J ct  ~  •  ^ 


^Tcr  =  •  1 


(3.42) 
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and  the  gain  would  be  stabilized  at  =  20.  Now,  if  a  large  input  occurs 

(pilot  command  input)  the  system  will  "jump"  to  the  locus  for  '^qct  ~  and 
with  the  gain  at  —  20  the  system  will  be  unstable  -0.3  for  the  actuator 

or  detected  mode).  To  avoid  this  instability  the  frequency  sensor  will  have  to 
decrease  very  rapidly  -  as  fast  as  the  large  input  increases  the  gain.  In 

a  "nut-sheli",  where  nonlinearities  in  the  control  system  cause  changes  in  the 
locus  of  the  mode  being  uced  to  adapt  the  system  (the  mode  whose  frequency 
is  the  input  to  the  frequency  sensor),  then  the  adaptive  loop  must  be  able  to 
respond  (change  the  gain)  fast  enough  to  accommodate  the  changes  in  allowable 
gain  which  result  from  the  locus  shift.  This  applies  even  if  the  shift  is  just 
one  of  gain  and  not  of  position  in  the  s-plane.  For  the  example  case,  the  adap¬ 
tive  loop  must  be  able  to  respond  at  least  as  fast  as  the  airplane  short  period 
and  perhaps  faster.  Conclusions  based  on  the  foregoing  linear  analysis  are 
verified  by  the  analysis  of  Section  4.  3  which  properly  includes  the  nonlinearities 
and  the  adaptive  loop  dynamics. 

The  damping  of  rate  gyros  tends  to  vary  with  temperature,  even 
for  "temperature  compensated"  rate  gyros.  A  typically  modest  variation 
for  such  a  rate  gyro  was  selected,  .2<  .  7 ,  and  its  effect  is  evaluated  in 

Figure  20.  The  root  locus  variations  are  extreme,  and  the  situation  presented 
in  Figure  20  is  clearly  unsuitable  for  application  of  the  G.E.  frequency  sensor. 
The  variation  of  the  gyro  poles  causes  far  too  extreme  a  variation  in  the  locus 
of  the  mode  that  goes  unstable  and  would  not  be  detected  by  the  frequency  sensor. 
Figure  20  is  a  good  illustration  of  one  of  the  basic  problems  involved  in  applying 
the  G.  E.  frequency  sensing  technique  to  adaptive  control.  The  root  locus  of  the 
detected  mode  must  be  well  defined  and  stabilized  (fixed  in  the  s-place).  This 
means  knowing  many  of  the  control  system  characteristics  quite  accurately,  es¬ 
pecially  those  with  natural  frequencies  near  the  frequency  to  be  detected  by  the 
frequency  sensor.  Furthermore,  the  characteristics  of  such  control  system 
elements  must  be  relatively  invariant.  If  elements  in  the  control  system  have 
varying  characteristics,  these  variations  must  be  considered,  just  as  are  the 
variations  in  airplane  characteristics. 
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3.4  Adaptive  Loop  Characteristics  j 

l  . 

Study  of  the  adaptive  loop  characteristics  requires  introducing  the  spe¬ 
cific  characteristics  of  the  frequency  sensor.  The  frequency  sensor  has  been 
simulated  on  both  the  EAS£  analog  computer  and  the  REAC  analog  computer.  A 
description  of  the  G.E.  frequency  sensor  is  given  in  the  following  paragraphs. 

A  functional  diagram  of  the  ana'.og  set-up  for  the  frequency  sensor  is 
given  in  Figure  9.  The  complete  EASE  analog  schematic  is  given  in  Figure  48, 
Appendix  C  and  the  REAC  analog  schematic  is  given  in  Figure  49  of  Appendix  D. 

These  analog  circuits  are  similar  to  that  used  by  G.E.  in  their  analog  analysis  of 
the  system.  The  frequency  sensor  operation  is  explained  using  the  functional 
diagram  Figure  9.  The  input  to  the  frequency  sensor,  St  ,  is  obtained  from  a 
signal  available  in  the  closed-loop  system  that  contains  the  frequency  1/  which 
is  to  be  stabilized  at  the  chosen  reference  frequency,  4-0  .  A  satisfactory  signal 
can  be  obtained  before  or  after  the  variable  gain  potentiometer,  fcj  in  Figure  8, 
or  from  the  output  of  the  canceller.  The  input  signal  S(  goes  through  a  relatively 
broad  bandpass  filter  hr.  J  where  X  is  chosen  so  that  the  peak  ampli¬ 
tude  of  the  filter  is  approximately  at  the  reference  frequency  ■£<>  .  This  filter 
produces  a  12  db/octave  rise  and  12  db/octave  decay  around  this  peak  fre¬ 
quency.  The  purpose  of  the  filter  is  to  isolate  the  frequency  of  the  mode  being 
detected  by  the  frequency  sensor,  from  either  lower  frequency  modes  such  as 
the  airplane  short  period,  or  from  higher  frequency  modes  such  as  structural 
modes  or  rate  gyro  modes.  The  signal  at  S ^  then,  is  largely  composed  of  fre¬ 
quencies  around  the  reference  frequency  •£#  .  On  the  block  diagram,  the  first 
half  period  of  the  signal  at  is  designated  as  being  equal  to  7} ,  where 

can  be  larger  or  smaller  than  £0 .  The  sinusoidal  signal  S2  is  converted 
to  »  square  wave.  The  zero  crossings  of  this  square  wave  signal  become,  the 
controlling  quantities  for  the  gain  changer  operation.  The  square  wave  is  now 
converted  to  a  pulse  signal  of  approximately  .  02  seconds  duration  at  the  time 
of  each  zero  crossing.  This  pulsed  signal  provides  the  control  inputs  for  the 
elements  in  the  lower  blocks  of  Figure  9. 

The  integrator,  at  the  lower  left  of  Figure  9,  has  an  adjustable  input, 
called  T0  adjust,  which  controls  the  rate  at  which  the  input  IN  sweeps  along  ; 

the  function  programmed  into  the  function  generator.  This  adjustment  sets  ;  - 

i 

L 
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the  frequency  fc  .  If  f^  is  present  in  3%  ,  no  change  is  produced  in  the  gain. 

Any  deviation  from  in  will  result  in  an  increase  cr  decrease  in  gain. 

The  integration  starts  approximately  at  the  end  of  each  gate  pulse  and  continues 
through  the  next  pulse.  The  integrator  is  reset  by  the  pulse  (integrator  reset 
signal)  derived  from  the  RC  network  and  diode  as  indicated  on  the  functional 
diagram.  The  output  of  the  function  generator  (OUT)  thus  depends  on  the  sweep 
rate  provided  by  the  integrator  and  the  time  each  pulse  occurs  in  the  integrator 
reset  signal.  The  gate  allows  only  the  last  .02  seconds  of  the  output  (OUT)  sig¬ 
nal  to  pass  (for  any  one  integration  cycle).  The  output  (OUT)  of  the  function  gen¬ 
erator  that  passes  through  the  gate  will  be  zero  if  the  pulse  occurs  at  the  set 
time  (T£  =  T  0  ),  positive  (+)  if  it  occurs  sooner  (T^  T0  ),  or  negative  (-) 

if  it  occurs  later.  The  specific  input  T£  shown  in  Figure  9  occurs  before  T0 
indicating  that  the  frequency  at  S2  is  too  high.  Small  deviations  from  the  ref¬ 
erence  frequency  produce  a  linear  output.  Large  deviations,  however,  are 
weighted  less  as  governed  by  the  function  generator. 

When  each  pulse  occurs,  the  gate  ie  opened  and  the  output  (OUT)  of  the 
function  generator  is  integrated  for  the  duration  of  the  pulse.  The  gain  control 
servo  positions  a  potentiometer,  a  variable  voltage  divider,  which  is  the  variable 
forward-loop  gain  of  the  closed-loop  system.  The  polarity  of  the  output  signal 
from  the  gate  determines  whether  the  gain  increases,  decreases,  or  remains 
constant.  The  maximum  increment  of  gain  change  that  can  occur  for  each  pulse 
duration  is  determined  by  the  peaxs  of  the  function  generator.  This  increment, 
which  ca.i  be  adjusted,  determines  the  speed  of  adaption.  Since  the  integrator 
only  operates  during  the  time  of  the  pulse,  the  servo  moves  in  steps.  These 
steps  can  be  filtered  by  slowing  down  the  response  of  the  servo;  however,  this 
filtering  also  will  decrease  the  speed  of  adaption.  These  aspects  of  the  fre¬ 
quency  sensor  are  discussed  further  is  Section  IV, 

Dynamic  responses  for  the  controlled  airplane  were  obtained  on  the  analog 
with  the  frequency  sensor  in  the  system.  Aircraft  (X- 15) .actuator  (second  order), 
and  control  system  characteristics  used  were  those  given  in  Section  3.3,  and  the 
flight  condition  studied  corresponded  to  t  =  90  sec  in  the  re-entry.  The  prefilter, 
for  obtaining  the  desirable  response,  was  included  in  the  input.  The  analog 
responses  obtained  are  presented  in  this  section,  along  with  brief  commentary. 

No  extensive  analysis  was  performed  on  the  data. 
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To  interpret  the  time  histories  of  the  system  gain,  it  is  neces¬ 
sary  to  explain  how  this  gain  was  mechanized  on  the  analog  computer.  Two 
stages  were  used.  The  first  stage  was  a  potentiometer,  or  attenuator,  controlled 
by  the  gain  changer  servo,  which  varied  the  input  to  a  fixed  gain  amplifier. 


The  variable  gain  can  be  described  by  the  following  equations  (Refer  to 
Figure  8): 


^3  = 


(3.43) 


(3.44) 


kl  ~  k!x  + 

X  "/> 


i  *r„s 


(3.45) 


*«“[**•  *775^]* 


(3.46) 


where  is  potentiometer  setting 

&  is  trimmed  or  initial  potentiometer  setting 
is  change  of  the  potentiometer  from  l&x0 
Kp  is  gain  of  amplifier 

£s  is  output  of  frequency  sensor  integrator 

is  time  constant.  This  can  be  thought  of  as  a 
filter  time  constant  for  ,  or  as  the  time 
constant  of  the  servo  which  moves  the 
potentiometer. 

From  these  equations,  it  is  seen  that  a  "unit"  change  in  the  output  of  the  fre¬ 
quency  sensor  integrator  (  Cg  )  causes  a  "unit"  change  in  the  potent- m  re  ter 
setting,  If  K*o  is  large,  that  is,  if  the  potentiometer  setting  is  near  1.0, 

a  "unit"  change  in  or  A**  will  cause  a  small  percentage  change  in  or 

.  However,  if  Kya  is  small,  the  "unit"  change  in  £>g  or  causes  a  large 
percentage  change  in  Al^.  or  (£g  . 
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The  analog  records  of  Figures  21  through  23  demonstrate  this  situation. 

Also,  these  analog  records  contain  data  on  the  effect  of  varying  the  time  constant, 

T’jj  .  The  analog  records  presented  in  Figure  21a  show  how  the  speed  of  adaption 

is  affected  by  Tx  when  the  gain  is  too  high  initially,  and  the  desired  gain  is 

such  that  K*  is  large.  The  same  effect  is  shown  in  Figure  21b  when  the  desired 
*0 

gain  is  such  that  K#  is  small.  Figure  21c  shows  the  case  when  the  desired 
gain  is  very  low,  near  zero.  Finally,  Figure  22  shows  both  the  case  of  initially 
too  high,  and  initially  too  low  .  In  all  of  these  figures,  step  command  inputs 

are  inserted  to  show  the  effect  of  adaptive  gain  dynamics  on  the  airplane  response. 
Also  shown  are  the  effects  of  pulses  forced  into  the  gain  changer  servo.  In 
general,  the  responses  show  good  characteristics  at  high  K*o  ,  and  poor  ones 
at  low  Kx0  .  At  low  %Xo  •  *ke  system  is  noisy  and  changes  ir.  £ ^  cause  severe 
disturbances  in  the  airplane's  response.  A  comparison  of  Figures  23a  and  23b 
shows  most  clearly  the  effect  of  changes  in  AKX  or  .  The  following  equation 
for  the  error  signals  in  the  system  helps  interpret  the  results  (see  Figure  8). 

<S.er  =  -  *>  ■#<*»„  **<x)  <3'47) 

A  change  in  Kx  ,  that  is,  a  ,  has  the  same  effect  as  the  change  in  £z 

(the  integral  of  the  error). 

The  inputs  from  &KX  could  be  made  more  tolerable  if  the  integration 
were  performed  after  the  gain  changer  rather  than  before.  In  this  case  would 
be  operating  on  the  error  (  £  )  and  not  the  integral  of  the  error  ( -  i-  6  ). 

The  error  is  usually  small,  but  the  integral  of  the  error  usually  has  a  large  value 
required  to  trim  the  aircraft.  Thus  the  "gain"  for  a  AKX  input  would  be  much 
smaller,  in  fact,  zero  except  during  transient  responses,  if  the  integration 
followed  the  gain  changer.  This  would  require  modification  of  the  method  of 
mechanizing  the  feedback  loops  of  the  control  system. 

Another  method  for  improving  the  performance  of  the  adaptive  loop 
would  be  to  mechanize  the  variable  gain,  Kj  »  such  that  a  "unit"  change  in  £3 
or  r  would  result  in  a  constant  percentage  change  in  the  existing  gain. 
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SECTION  IV 

STUDY  OF  SPECIFIC  PROBLEMS  IN  APPLICATION  OF 

G.E.  SYSTEM 


4.  1  Introduction 

The  analysis  of  the  G.E.  System  discussed  in  Section  III  rjed  linear 
analysis  techniques  based  principally  on  the  root  locus  method,  af>d  simula¬ 
tion  of  the  system  on  an  analog  computer.  The  X-15  was  the  study  vehicle, 
and  the  equations  and  parameters  used  are  described  in  Appendix  A. 

The  problems  examined  in  the  analysis  of  Section  III  may  be  listed  as 
follows : 

1.  Model  concepts 

2.  Selection  of  system  parameters 

3.  Responses  to  step  pitch  rate  commands 

4.  Response  to  gusts 

5.  Modification  of  the  basic  airplane  characteristics  - 
particularly  unstable  configurations 

6.  Sensor  and  actuator  dynamics 

7.  Dynamics  of  the  adaptive  loop. 

The  analog  computer  used  in  the  above  studies  was  a  small  EASE  com¬ 
puter,  and  simulation  of  the  airplane,  the  actuator,  the  reference  model  and 
the  frequency  sensor  used  it  to  full  capacity.  .Few  components  were  available 
for  simulation  of  any  additional  elements  or  nonlinearities.  Since  the  results 
of  the  root  locus  studies  had  indicated  that  actuator  non-linearities  could  be  a 
very  important  problem  area  for  the  G.E.  System,  it  was  decided  to  remech¬ 
anize  the  problem  on  a  larger  REAC  computer  so  that  actuator  nonlinearities 
could  be  studied.  This  was  done  and  the  following  problems  were  examined 
using  this  machine. 

1.  The  response  of  the  system  to  inputs  at  the  actuator  valve 

2.  The  effect  of  actuator  nonlinearities 
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3.  Characteristics  of  the  frequency  sensor 

4.  The  effect  of  "noise"  in  the  rate  gyro  output 


The  analog  schematics  of  the  problem  as  it  was  mechanized  on  the  REAC 
analog  computer  are  contained  in  Appendix  D. 

4.2  Response  to  Pilot  Inputs  at  the  Actuator  Valve 

The  method,  illustrated  in  Figure  52,  by  which  the  valve  servo  displace¬ 
ments  and  pilot  commands  from  the  mechanical  system  can  be  mixed  and  applied 
to  the  power  actuator  valve  in  the  X-15,  suggested  that  it  might  be  of  interest  to 
determine  the  response  of  the  system  to  pilot  inputs  through  the  mechanical  sys¬ 
tem  to  the  power  actuator  valve. 


The  block  diagram  of  Figure  24  illustrates  the  situation  in  a  manner  suit¬ 
able  for  analysis.  The  elevator  command  signal  can  be  moved  to  the  left  of  the 
summing  point  by  multiplying  by  the  inverse  of  the  transfer  functions  of  the  boxes 
between  the  summing  point  and  the  power  actuator.  From  this  diagram,  it  can 
be  observed  that  a  step  command  will  now  result  in  a  step  pitch  angle  response 
rather  than  a  step  pitch  rate  response.  Also,  the  response  amplitude  will  now 
be  dependent  on  the  value  of  ,  i.  e. ,  a  larger  value  will  result  in  a  smaller 
pitch  response  for  a  given  elevator  command.  In  equation  form: 
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Q  Te  S+t 


ignoring  the  high  frequency  zeros  due  to  C.he  valve  servo. 

Similarly,  the  angle  of  attack  response  can  be  written: 


(4.3) 
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The  angle  of  attack  response  to  an  elevator  command  will  consist  of  a 
fast  rise  and  then  a  decay  to  zero  at  a  rate  dependent  cn  the  value  of  .  How¬ 
ever,  for  the  flight  conditions  where  the  dynamic  pressure  is  very  low,  the  value 
of  is  very  large.  This  in  effects  puts  the  pole  (W  at  the  origin  and  can¬ 
cels  the  zero  due  to  the  (S)  in  the  numerator.  Under  these  conditions,  would 
be  very  small  and  could  be  ignored  along  with  the  zeros  due  to  the  valve  servo. 
The  response  would  be 
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Thus  one  would  expect  a  steady-state  angle  of  attack  in  response  to  a  step  ele¬ 
vator  command  for  the  high  altitude  low  ^  flight  conditions. 

The  analog  records  of  Figure  35  support  the  above  conclusions.  Re¬ 
sponses  of  the  open-loop  airplane  to  elevator  steps  are  shown  in  Figure  25a  for 
flight  conditions  of  t  =  0,  40,  and  90  seconds  in  the  re-entry.  The  records  of 
Figure  25b  are  responses  of  the  X-15  with  the  G.S.  system  to  pitch  rate 
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commands  for  flight  conditions  of  t  =  40  and  90  seconds,  and  finally,  in  Figure 
25c  are  responses  with  the  G.E.  System  to  elevator  commands  for  flight  condi¬ 
tions  of  t  =  0,  40,  and  90  seconds.  The  effect  of  the  value  of  on  the  response 
to  elevator  commands  can  be  seen  in  Figure  25c  by  comparing  the  responses  for 
t  =  90  with  those  for  t  =  0  and  40.  The  effect  of  Ky  on  the  amplitude  can  be  seen 
by  comparing  tv-o  responses  at  the  same  flight  condition.  Unfortunately,  the 
same  recording  scales  were  not  used  in  all  of  the  above  records,  and  this  tends 
to  mask  the  effect  of  My  on  the  response  amplitude  in  some  of  the  records. 

The  two  responses  in  Figure  25b  to  pitch  rate  commands  at  t  -  90  flight 
condition  for  two  gain  values,  My  =  2  and  4,  illustrate  the  effect  on  the  O  re¬ 
sponse  of  not  having  the  gain  high  enough  to  drive  the  integrator  pole  to  the  air¬ 
plane  zero  at  as  was  discussed  in  Section  III,  the  effect  being  that  the  re¬ 

sponse  is  more  like  a  third-order  system  than  the  second-order  model.  The  re¬ 
sponse  to  pitch  rate  commands  for  the  flight  condition  of  t  =  40  in  Figure  25b 
illustrates  the  extremely  slow  angle  of  attack  response  obtained  when  is 

large.  The  pulse  inputs  to  the  forward  loop  in  Figure  25b  illustrate  the  damping 
of  the  actuator  mode  for  the  two  gain  levels. 

In  summary,  it  was  found  that  the  effect  of  applying  pilot  commands 
through  the  mechanical  stick  to  the  elevator  actuator  valve  rather  than  as  pitch 
rate  commands  was  to  change  the  system  from  a  pitch  rate  command  system  to 
a  pitch  attitude  command  system  with  the  amplitude  of  the  response  being  an 
inverse  function  of  the  My  gain. 

After  the  analog  study  reported  here  had  been  completed.  References  22 
and  23  were  obtained.  It  is  interesting  to  note  that,  in  the  X-15  simulator  study 
of  the  G.E.  Syste^n  reported  in  these  references,  the  pilot  input  to  the  system 
was  a  combination  of  th.  normal  mechanical  input  to  the  actuator  valve  and  an 
electrical  pitch  rate  command.  This  input  can  oe  expressed  as  follows: 
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Thus,  by  moving  the  stick,  the  pilot  commanded  a  combination  of  pitch  rate 
and  pitch  attitude  with  the  pitch  attitude  command  inversely  proportional  to  the 
gain. 

4.  3  Effect  of  Actuator  Nonlinearities 
4.  3.  1  Root  Locus  Studies 

The  root  locus  studies  of  Section  III  indicated  that  a  nonlinearity 
which  caused  the  actuator  time  constant  to  change  with  amplitude  v/ould  result 
in  changes  in  the  root  locus  of  the  actuator  mode  that  could  have  serious  conse¬ 
quences  on  the  operation  of  the  G.  E.  frequency  sensor  and  gain  changer. 

Further  root  locus  studies  were  made  using  the  ESIAC  root  locus 
plotter  to  establish  gain  levels  and  frequencies  of  the  actuator  mode  for  various 
flight  conditions  and  two  actuator  time  constant,  TAcr  =  .  5  and  .  1  seconds.  The 
results  of  these  studies  are  illustrated  in  Figure  26a  through  d.  It  was  noted 
for  a  flight  condition  of  t  =  74  and  the  slow  actuator  that  the  root  locus  plotter 
was  undecided  as  to  the  locus  path  of  the  airplane  poles  and  the  actuator  locus 
at  certain  gains.  The  frequency  of  the  airplane  short-period  poles  was  increased 
and  more  root  locus  plots  were  made.  It  was  found  that  a  locus  shift  would 
occur  with  the  result  that  the  aircraft  poles  would  go  unstable  and  the  actuator 
locus  would  go  to  the  model  zeros  as  shown  in  Figure  26d.  However,  there 
was  no  significant  change  in  the  frequency  at  which  the  loop  became  unstable. 

The  values  of  /Cj  gain  ard  actuator  mode  frequency  for  =0  and  .  2  obtained 
from  the  ESIAC  studies  are  tabulated  in  Table  IV.  Similar  data  for  g  -  0 
were  obtained  from  the  analog  computer  by  pulsing  the  forward  loop  at  various 
values  of  J¥j  until  the  value  of  iCj  was  reached  where  the  transient  response 
had  zero  damping.  The  frequency  of  the  oscillation  could  be  read  from  the 
time  history.  This  data  is  tabulated  in  Table  V.  The  agreement  is  fairly  good 
for  flight  conditions  of  t  =  60,  74,  and  9C  seconds,  but  for  the  low  dynamic  pres¬ 
sure  conditions  of  t  =  0,  20,  and  40  seconds,  the  agreement  with  the  ESIAC 
results  is  not  so  good.  It  is  thought  that  this  is  due  to  inaccuracies  in  the 
analog  solutions.  The  problem  was  scaled  for  the  high  dynamic  pressure  flight 
conditions  and  at  the  low  fa  conditions,  the  analog  voltages  were  too  low  for 
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accurate  solutions.  It  is  of  interest  to  note  that  during  the  checkout  of  the  analog 
set-up,  in  two  instances,  changes  in  the  mechanization  (to  pick  up  gain  or  addition 
of  a  small  capacitor  in  an  amplifier  feedbackto  reduce  noise)  introduced  unwanted 
poles  in  the  frequency  band  of  interest.  The  result  was  that  the  root  loci  were  shifted 
and  the  gain  and  frequency  limits  for  stability  were  changed.  This  is  reported 
as  an  example  of  the  difficulties  that  practical  considerations  can  introduce.  It 
emphasizes  the  point  that  in  applying  the  G.  E.  System  to  a  control  problem,  it 
is  necessary  that  all  poles  and  zeros  in  the  frequency  band  that  might  affect  the 
root  loci  be  known. 

The  gain  (ESLAC  data)  required  to  maintain  the  actuator  mode  damping 
at  £  =  .2  for  each  of  the  two  actuator  time  constants  is  plotted  in  Figure  27. 

4.3.2  Analog  Simulation  of  Components 

As  part  of  the  analog  check-out,  and  to  obtain  accurate  information 
for  the  ESIAC  studies,  frequency  response  measurements  were  made  on  the  net¬ 
works  used  to  represent  the  rate  gyro  and  valve  servo  together  with  the  G.  E. 
feedback  network.  The  results  of  these  tests  are  illustiated  in  Figures  28  and 
29.  The  rate  gyro  was  represented  by  a  second-order  system  with  a  natural 
frequency  of  13.2  cps  and  a  damping  ratio  of  £  =  .68.  The  valve  servo  was 
represented  by  a  second-order  system  with  a  natural  frequency  of  35  cps  and  a 
damping  ratio  of  £  =  .43.  The  frequency  of  this  servo  is  probably  higher  than 
the  servo  actually  installed  in  the  X-15,  but  this  is  not  thought  to  be  very  criti¬ 
cal  to  the  studies  made.  The  main  reason  for  using  these  servo  dynamics  was 
that  the  pad  had  been  built  up  for  another  program  and  was  available. 

4C3.3  Nonlinear  Actuator 

Information  available  concerning  the  X-15  longitudinal  control  sys¬ 
tem  indicated  that  its  response  rate  was  a  function  of  amplitude  and  frequency, 
i.e.,  act  *  *5  seconds  for  low  amplitude  and  T4ct  *  .  1  seconds  for 

large  amplitudes.  This  nonlinearity  was  considered  to  be  primarily  due  to  the 
flow  characteristics  of  the  actuator  control  valve.  Additional  nonlinearities 
that  could  be  included  were  actuator  rate  and  displacement  limits  and  slop  or 
backlash  in  the  control  linkages.  In  the  studies,  all  of  these  nonlinearitiea 
except  backlash  were  simulated. 
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A  curve  of  rate  of  change  of  flow  with  valve  displacement  which  con¬ 
tained  the  nonlinear  characteristics  to  be  simulated  was  found  in  Reference  18 
and  is  considered  to  be  representative.  This  curve  is  replotted  in  Figure  30. 

The  scales  were  normalized  and  the  curve  was  reshaped  at  the  bottom  to  give  a 
maximum  value  of  =  10  and  minimum  value  of  dQ/a^  =  2.  These  values 

of  correspond  to  actuator  time  constant  of  T?aot  -  ♦  1  and  .  5,  respectively. 


The  curve  of  in  Figure  30  was  numerically  integrated  to  obtain 

data  which  was  replotted  as  flow,  Q  ,  vs.  valve  displacement,  /y  ,  in  Figure 
31.  A  function  generator  was  built  to  produce  this  relation  as  indicated  in  the 
block  diagram  of  the  nonlinear  actuator,  Figure  32.  By  proper  adjustment  of 
the  limiters  on  the  output  of  amplifier  IB,  limits  are  imposed  on  the  actuaior 
rate  and  by  adjustment  of  the  limiters  on  the  output  of  the  integrator  5B,  limits 
are  imposed  on  the  actuator  displacement.  A  detailed  circuit  schematic  is  con¬ 
tained  in  Appendix  D.  The  position  limits  used  were  § e  =  15°  and  -35°,  the 
rate  limits  were  varied  from  a  nominal  value  of  =  20° /sec.  The  frequen¬ 

cy  response  and  transient  responses  of  Figures  33  and  34  illustrate  the  non¬ 
linear  response  of  the  power  actuator  with  valve  displacement  amplitude  ob¬ 
tained  in  the  analog  simulation. 


4.3.4  Analog  Records  of  X-15  with  G.E.  System 
Including  Nonlinear  Actuator 

The  purpose  of  the  analog  studies  was  to  explore  and  to  demonstrate 
the  posrible  effects  of  nonlinearities  on  the  operation  of  the  G.E.  frequency  sen¬ 
sor  and  gain  changer.  From  the  linear  analysis  based  on  the  root  locus  techniques, 
it  was  found  that  the  possibility  existed  that  during  conditions  when  small  actuator 
valve  displacements  were  required,  the  gain  changer  would  select  a  gain  level 
consistent  with  the  root  locus  for  the  slow  actuator.  If  then  a  large  maneuvering 
command  were  inserted  or  a  gust  encountered,  the  system  would  "jump"  to  the 
locus  for  the  fast  actuator  and  the  gain  would  be  too  high  by  a  factor  of  about  2.5, 
and  the  system  would  go  unstable.  Since  the  gain  would  be  too  high  by  a  consid¬ 
erable  factor,  the  damping  ratio  of  the  unstable  mode  would  be  a  large  negative 
value  and  the  instability  would  be  expected  to  be  nearly  explosive.  However,  the 
actuator  rate  limits  prevent  this  from  happening  and  the  result  is  a  fixed  amplitude 
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limit  cyclle,  the  amplitude  of  which  is  a  function  of  the  system  loop  gain.  Thus, 
the  rate  limit  has  a  beneficial  effect  in  that  it  restricts  the  amplitude  of  unstable 
oscillations.  Another  phenomenon  characteristic  of  limit  cycles  which  has  spe¬ 
cial  significance  for  the  operation  of  the  G.E.  frequency  sensor  is  that  the  fre¬ 
quency  of  the  limit  cycle  decreases  with  increasing  gain.  Thus,  the  possibility 
exists  that  the  valve  flow  nonlinearity  could  cause  the  gain  to  be  driven  too 
high  resulting  in  an  unstable  oscillation  which  would  be  converted  into  a  limited 
amplitude  limit  cycle  at  a  frequency  lower  than  the  free  oscillation.  If  the  fre¬ 
quency  of  this  limit  cycle  falls  below  the  reference  frequency  of  the  G.E,  fre¬ 
quency  sensor,  the  gain  changer  would  drive  the  gain  to  even  higher  levels  and 
a  divergent  situation  would  result.  The  above  described  operation  has  been 
demonstrated  on  the  analog  computer  and  is  illustrated  by  the  sample  reccrus 
reproduced  in  Figures  35,  36a,  and  36b. 

In  Figure  35,  the  operation  of  the  system  with  a  linear  actuator  (%lers 
0. 1  sec)  at  flight  condition  t  =  90  is  illustrated.  In  these  records,  the  input  to 
the  frequency  sensor  was  taken  from  the  output  of  the  "canceller"  or  washout, 
portion  of  the  G.E.  feedback.  In  the  absence  of  "pilot"  inputs  to  the  elevator, 
the  frequency  sensor  had  a  tendency  to  drive  the  gain  up.  This  could  have  been 
the  result  of  low  frequency  noise  or  the  integrator  of  the  frequency  sensor  may 
have  had  a  bias.  In  any  event,  the  gain  would  increase  until  a  disturbance  ex¬ 
cited  the  actuator  mode  and  then  the  frequency  sensor  would  reduce  the  gain  and 
the  process  would  tend  to  repeat.  The  amplitude  of  the  oscillations  was  small 
(  Sf  <  ±.25°)  and  the  rate  limits  were  not  exceeded.  If  continuous  small  ampli¬ 
tude  elevator  disturbances  were  applied  by  the  "pilot",  the  system  would  tend  to 
maintain  a  more  constant  gain  level.  For  the  record  of  Figure  35,  the  reference 
frequency  was  set  at  0)0  =  5  cps. 

In  Figure  36a  the  effects  of  the  nonlinear  actuator  valve,  the  actuator 
rate  limits,  and  the  setting  of  the  reference  frequency  of  the  frequency  sensor 
are  demonstrated.  At  the  beginning  of  the  record,  in  the  absence  of  gusts  or 
pilot  inputs,  the  gain  changer  increased  the  gain  as  permitted  by  the  slow  actu¬ 
ator  time  constant,  '**ct  =  -5.  which  is  stable  for  gains  of  A'j  =  18.  When 
the  gain  had  increased  to  tfy  ts  17,  however,  the  system  became  disturbed  and 
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developed  valve  displacements  large  enough  to  decrease  the  actuator  time  con¬ 
stant  to  near  t -  .  1  for  which  the  maximum  stable  gain  is  only  =  4.4. 

Tjus  the  system  became  suddenly  very  unstable  and  went  into  a  limit  cycle, 
which  in  this  case  had  a  frequency  slightly  higher  than  the  reference  frequency 
of  <u0  st  S  i  cps  and  the  gain  was  driven  down.  Close  inspection  of  this  record 
indicates  a  change  in  the  amplitude  of  the  signal  from  the  frequency  sensor  gate 
amplifiers,  and  a  corresponding  change  in  the  rate  at  which  the  gain  decreased, 
initially  as  the  free  oscillation  was  developing,  it  had  a  considerably  higher  fre¬ 
quency  than  the  reference  frequency,  and  the  gate  opened  when  the  output  of  the 
function  generator  had  nearly  its  maximum  value.  As  the  limit  cycle  developed, 
the  frequency  of  the  oscillation  decreased  and  the  gate  then  opened  when  the  out¬ 
put  of  the  function  generator  was  lower  or  near  the  zero  crossing  and  the  gain 
changer  slowed  down.  The  result  was  that  it  took  more  tha.t  ten  seconds  for  the 
system  to  recover.  The  recorder  was  then  stopped  and  the  reference  frequency 
was  increased  to  approximately  c00  -  5.2  cps.  Again,  the  gain  was  increased 

by  the  frequency  sensor  until  a  disturbance  set  off  the  instability.  However,  this 
time  the  frequency  of  the  fully  developed  limit  cycle  was  lower  than  the  reference 
frequency  and  the  gain  changer  increased  the  gain  to  the  maximum  permitted  by 
the  servo.  The  activ  tor  rate  limits  were  set  to  approximately  £«  =  ±20* /sec 
f^r  these  records.  The  unstable  gain  changer  operation  could  be  made  to  occur 
with  -  5  cps  by  reducing  the  rate  limit  to  about  ±10°/sec,  but  the  most  cri¬ 
tical  parameter  i?  the  frequency  sensor  reference  frequency  setting. 

Figure  2Sb  illustrates  the  operation  of  the  system  fora  flight  condition  of 
t  =  40.  For  this  record,  the  maximum  for  th-2  slow  actuator  was  determined 
and  the-  feedback  potentiometer  of  the  # 5  amplifier  was  adjusted  to  permit  this 
gain  (  -  175)  whan  the  servo  driven  attenuator  was  at  the  high  end.  The  ac¬ 

tuator  rate  'imits  were  set  at  app .-oximately  ±25*/sec  and  the  reference  fre¬ 
quency  was  set  at  <r>0  =5.3  cps.  As  before,  the  gain  increased  until  a  dis¬ 
turbance  caused  the  oscillation  to  develop.  The  first  two  times,  a  disturbance 
occurred  when  the  gain  was  not  too  high  and  the  system  recovered.  However, 
the  diird  time  the  disturbance  did  net  occur  until  the  gain  had  exceeded  the  cri¬ 
tical  value  that  caused  the  limit  cycle  frequency  to  be  lower  than  the  reference 
frequency  and  the  frequency  sensor  caused  the  gain  to  increase  to  the  limit. 


58 


ASD-TR-61-104,  Vol.  I 


In  summary,  it  has  been  demonstrated  that  although  the  G.  E.  System 
with  its  frequency  sensor  and  gain  changer  is  capable  of  keeping  the  closed-loop 
performance  nearly  constant,  independent  of  airplane  characteristics  and  flight 
condition,  the  system  can  encounter  real  difficulties  if  the  actuator  has  non- 
linearities  or  changes  its  operation  in  such  a  manner  as  to  cause  the  root  locus 
of  the  node  which  becomes  unstable  to  shift.  In  the  case  of  the  system  studied 
here,  the  critical  parameters  are  the  actuator  time  constant,  "t^cr  *  the  can¬ 
celler  time  constant,  Tq  ,  and  the  rate  gyro  dynamics.  A  change  in  any  of 
these  parameters  with  environment  or  time  can  cause  changes  in  the  root  loci 
such  that  the  frequency  sensor  and  gain  changer  may  not  operate  properly.  Thus, 
as  pointed  out  in  Reference  20,  the  adaptive  control  system  does  not  require 
specific  information  on  the  exact  configuration  cr  flight  condition  other  than  the 
over-ali  range  of  the  variables.  However,  the  behavior  of  the  system  is  often 
quite  strongly  dependent  on  the  actuator  characteristics.  This  situation  has 
the  advantage  that  the  required  system  adjustment  can  probably  be  determined 
on  the  ground  prior  to  the  flight. 

4.4  Frequency  Sensor  Characteristics 

Some  insight  into  the  characteristics  of  the  G.  E.  frequency  sensor  has 
been  acquired  through  the  process  of  mechanization,  check-out  and  operation 
of  the  device  on  the  analog  computer.  A  description  of  the  G.  E.  frequency  sen¬ 
sor  anc.  the  essence  of  how  it  works  was  included  in  Section  III.  Figure  9  is  a 
functional  block  diagram  of  the  frequency  sensor.  The  REAC  analog  schematic 
is  contained  in  Appendix  D.  Certain  characteristics  of  the  system  are  dis¬ 
cussed  in  the  following  paragraphs. 

First,  the  characteristics  of  the  bandpass  filter  will  be  considered. 

The  purpose  of  this  filter  is  to  isolate  the  frequency  component  of  the  input  sig¬ 
nal  that  is  due  to  the  mode  being  detected  by  the  frequency  sensor  from  lower 
frequency  components  due  to  pilot  inputs  or  airplane  short-period  oscillations 
and  higher  frequency  components  that  may  be  due  to  structural  modes  or  rate 
gyro  noise.  In  practice,  the  filter  does  indeed  function  in  this  manner  provided 
there  is  a  component  of  the  input  signal  that  is  due  to  the  actuator  mode.  However, 
if  the  gain  is  low  and  the  actuator  mode  has  relatively  high  damping  or  if  th“re 
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are  no  inputs  to  the  control  system  to  excite  the  actuator  mode,  then  there 
will  not  be  any  component  in  the  frequency  sensor  input  signal  due  to  the  ac¬ 
tuator  mode.  In  this  situation}  the  frequency  sensor  will  be  triggexcd  by  any 
noise  that  causes  a  aero  crossing  of  in  Figure  9*  Obviously,  zero  cros¬ 
sings  of  will  depend  on  the  relative  amplitudes,  weighted  by  the  bandpass 
filter,  of  the  various,  harmonic  components  of  the  input  signal  $,  .  Data  were 
taken  during  the  analog  computer  studies  to  illustrate  this  characteristic  and 
the  results  are  discussed  in  Section  4,  5. 

The  frequency  sensor  gates  are  opened  by  zero  crossings  of  the  signal 
S2  and  remain  open  for  a  fixed  length  of  time,  and  then  close.  The  length  of 
time  that  the  frequency  sensor  gate  is  held  open  is  somewhat  arbitrary  and  so, 
for  most  of  this  study,  the  gate-open  time  was  taken  as  .02  seconds  because 
this  was  the  value  G.  E.  used  in  their  flight  tests  and  analog  studies.  It  is  in¬ 
teresting  to  note  that  in  the  drawings  used  to  illustrate  the  function  of  the  sys¬ 
tem  in  Reference  20,  the  gate  is  drawn  as  about  5%  of  the  half  period  of  the 
reference  frequency,  whereas  actually,  .  02  seconds  is  nearer  to  20%  of  the 
half  period  ox  the  reference  frequency.  Thus,  for  u>a  #  30  rad/sec,  the  gate 

is  open  about  20%  of  the  time.  The  limits  on  gate-open  time  are  related  to  the 
shape  of  the  function  generator  curve  and  the  precision  with  which  the  gates 
can  be  made  to  open  and  close.  The  frequency  sensor  was  first  mechanized 
with  a  gate-open  time  of  about  5%  of  the  reference  frequency  half  period,  but  it 
was  found  that  this  resulted  in  unsymmetrical  outputs  from  the  two  gates  due  to 
tolerance  on  the  gate  width  and  any  small  bias  due  to  diode  characteristics. 

Thus,  it  appears  better  to  use  a  wide  gate  to  minimize  these  effects.  The  max¬ 
imum  gate  width  is  dependent  to  some  extent  on  the  shape  of  the  function  gen¬ 
erator  curve.  This  is  illustrated  in  Figure  37.  It  is  obvious  from  this  figure 
that  if  the  gate  is  held  open  toe  long,  then  the  significance  of  the  shape  of  the 
function  generator  is  lost.  That  is,  presumably,  the  function  generator  is 
shaped  in  a  particular  manner  so  as  to  weight  different  frequency  errors  dif- 
ferenly  and  if  the  gate-open  time  is  made  too  long,  the  effect  of  the  function 
shape  will  be  compromised.  Also  in  the  limit,  the  half  period  of  the  input  fre¬ 
quency  will  become  leas  than  the  gate  width  and  the  gates  will  be  open  continuously. 
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The  gate -open  time  also  has  an  effect  on  the  rate  at  which  the  gain 
changer  servo  runs.  This  is  because  the  gate  outputs  are  roughly  trains  of 
pulses  and  obviously,  if  the  pulses  are  made  wider,  the  integral  of  the  pulse 
train  for  a  given  period  of  time  will  be  larger.  The  servo  rate  can,  of  course, 
be  controlled  by  other  gains  in  the  system  so  the  effect  of  pulse  width  on  gain 
servo  rate  is,  for  the  most  part,  incidental. 

The  role  of  the  frequency  sensor  function  generator  is  to  weight  the 
frequency  errors  such  that  small  errors  ±lcps  from  the  reference  frequen¬ 
cy  will  result  in  large  error  signals  to  the  gain  changer  servo  while  very  large 
frequency  errors  will  be  weighted  less  by  the  function  generator  to  prevent 
random  noise  inputs  from  driving  the  gain  servo.  For  random  noise  inputs, 
the  result  is  as  described,  but  if  the  noise  is  not  random,  but  is  due  to  a  per¬ 
iodic  source  such  as  a  vibration  mode  or  electrical  pic’* -up,  then  the  frequency 
sensor  output,  i.e.,  the  rate  of  change  of  the  gain  servo,  becomes  not  only  a 
function  of  the  shape  of  the  function  generator  curve,  but  also  a  strong  function 
of  the  frequency  of  the  input.  This  situation  is  illustrated  in  Figures  37  and  38. 
In  Figure  37,  the  pulse  trains  are  drawn  for  input  frequencies  of  £c>/  =8.6  cps 
and  0*i  =3.5  cps.  In  both  cases,  the  gate  opens  on  the  function  generator  tails 
and  so  the  pulse  heights  and  widths  are  the  same.  However,  at  Ci>i  =8.6  cps 
the  pulses  are  much  closer  together  than  they  are  at  a)/  =3.5  cps,  and  so  the 
integral  of  these  pulses  will  drive  the  gain  changer  servo  at  a  higher  rate. 

Figure  38  is  a  plot  of  the  time  required  for  the  gain  servo  to  run  full  scale  in 
response  to  various  input  frequencies.  This  plot  shows  the  effect  of  input  fre¬ 
quency,  the  effect  of  the  function  generator  shape,  and  the  effect  cf  having  the 
input  half  period  equal  to  or  less  than  die  gate  width.  This  data  was  obtained 
from  analog  records  similar  to  those  of  Figure  39.  The  analog  traces  of  Figure 
39  also  illustrate  the  effect  of  input  frequency  on  the  gain  servo  rate.  Two  sets 
of  traces  are  included  in  Figure  39,  one  set  where  the  gate  outputs  are  summed 
and  then  integrated  and  another  set  where  the  gate  outputs  are  summed  and 
lagged  by  a  filter  with  a  one-second  time  constant  and  then  integrated.  In  both 
cases,  the  gate  outputs  are  in  steady-state  when  the  integrator  is  turned  on  30 
that  the  gain  servo  rate  obtained  is  the  steady  state  rate  and  does  not  contain  the 
transient  that  would  be  caused  by  the  lag  network.  This  effect  will  be  illustrated 
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by  records  in  Section  4.5.  The  records  of  Figure  39  illustrate  the  following: 
the  effect  of  the  input  frequency  on  the  servo  steady-state  rate,  the  effect  of 
the  function  generator  peaks  on  the  servo  steady-state  rate,  the  effect  of  the 
lag  network  in  smoothing  the  gate  outputs,  and  the  fact  that  the  lag  network  has 
no  effect  on  the  steady-state  rate  of  the  gain  servo. 

The  lag  network  was  placed  ahead  of  the  integrator  to  permit  easier 
control  of  the  servo  rate  and  to  keep  voltage  levels  within  the  limits  of  the  amp¬ 
lifiers.  See  Figure  49  of  Appendix  D. 

The  maximum  rate  of  change  of  the  gain  servo  for  small  frequency 
errors,  i.e. ,  on  one  peak  of  a  function  generator,  was  set  up  to  be  approxi¬ 
mately  40* /sec,  where  this  rate  refers  to  degrees  of  potentiometer  travel  on  a 
350*  potentiometer.  Thus,  the  servo  could  run  from  maximum  gain  to  mini¬ 
mum  gain  in  approximately  8-9  seconds. 

Another  observation  of  interest  is  that  a  bias  or  amplifier  drift  intro¬ 
duced  after  the  gates  will  result  in  a  tendency  for  the  gain  servo  to  drive  in  one 
direction.  If  die  bias  drives  the  servo  toward  high  gain,  the  actuator  mode  will 
become  lightly  damped  and  will  become  excited,  producing  an  input  to  the  fre¬ 
quency  sensor  which  will  cause  the  gain  to  be  corrected.  This  effect  is  illus¬ 
trated  in  Figure  35.  If  the  bias  or  drift  is  toward  low  gain,  the  system  charac¬ 
teristics  approach  those  of  the  open-loop  airplane  and  there  may  be  no  compo¬ 
nent  of  the  input  to  die  frequency  sensor  that  is  near  o>0  .  In  tkis  situation,  the 
gain  changer  will  operate  on  frequencies  outside  the  center  band  of  the  bandpass 
filter  and  the  center  portion  of  the  function  generator.  The  result  may  well  be 
that  the  frequency  sensor  will  not  recover  the  gain  to  the  proper  value  and  the 
system  response  will  not  be  that  of  the  model  at  all. 

In  summary,  die  G.  E.  frequency  sensor  operates  on  zero  crossings  of 
the  input  signal  after  its  various  harmonics  have  been  weighted  by  the  bandpass 
filter.  These  zero  crossings  open  the  gates  and  permit  signal  pulses  to  be  emitted 
which  have  their  amplitude  weighted  by  the  function  generator.  The  gain  changer 
servo  is  driven  by  a  signal  proportional  to  the  integral  of  the  pulse  train  from 
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the  frequency  sensor  gates  and  its  rate  is  therefore  a  function  not  only  of  the 
pulse  amplitudes,  but  also  of  their  frequency  of  occurrence.  These  charac¬ 
teristics  are  of  significance  in  determining  the  response  of  the  frequency  sen¬ 
sor  to  periodic  noise  in  its  input  signal. 

4.  5  Effect  of  Noise  in  the  Rate  Gyro  Output 

The  characteristics  of  the  G.E.  frequency  sensor  described  in  the  pre¬ 
ceding  section  suggest  the  possibility  that  periodic  noise  in  the  input  signal  to 
the  frequency  sensor  might  result  in  unsatisfactory  operation  of  the  gain  changer. 

To  investigate  the  e.'fects  of  periodic  noise  on  the  operation  of  the  fre¬ 
quency  sensor,  it  was  assumed  that  the  rate  gyro  was  mounted  such  that  it  could 
respond  to  structural  vibrations,  engine  vibration  or  any  other  mechanical  or 
electrical  oscillation  with  a  frequency  in  the  band  from  zero  to  30  cps.  Analog 
records  were  taken  of  the  system  operation  with  an  oscillatory  signal  injected 
into  the  output  of  the  amplifier  used  to  simulate  the  rate  gyro.  The  input  to  the 
frequency  sensor  was  taken  from  the  washout  network  or  canceller  part  of  die 
G.E.  feedback  system.  The  amplitude  of  the  oscillatory  input  was  adjusted  to 
equal  i. 05* /sec  of  pitch  rate  and  records  were  taken  at  various  frequencies 
from  3.5  to  30  cps.  This  data  was  used  to  obtain  the  plot  of  Figure  38.  The 
sample  analog  record  of  Figure  40a,  for  flight  condition  t  ~  90  and  the  linear 
actuator  with  =0.1  sec,  illustrates  the  operation  of  the  system  when  the 

noije  frequency  was  lower  than  the  reference  frequency.  The  gain  servo  was 
driven  to  high  gain  by  the  noise.  The  actuator  mode  became  light  1”  damped 
and  started  to  oscillate,  causing  an  input  to  the  frequency  sensor  with  a  frequen¬ 
cy  in  the  center  range  of  the  bandpass  filter.  The  bandpass  filter  amplified  the 
actuator  mode  frequency  and  attenuated  the  noise  input  so  that  the  zero  crossings 
of  $£  on  Figure  9  were  due  to  the  actuator  mode  and  the  function  generator  out¬ 
put  was  such  that  the  gain  was  reduced.  As  the  gain  decreased,  the  damping  of  the 
actuator  mode  increased  and  the  oscillation  due  to  the  actuator  mode  decayed. 

The  rate  gyro  noise  then  again  became  predominant  and  the  gain  was  increased 
to  start  another  cycle.  The  period  of  this  cycle  of  events  was  about  11  seconds. 

Figure  40  illustrates  the  operation  of  the  system  when  the  noise  frequency 
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was  higher  than  the  reference  frequency.  In  this  record,  the  noise  amplitude 
was  reduced  *by  a  factor  of  10  to  *.005*/sec  of  pitch  rate  at  8.5  cps.  The  amp¬ 
litude  of  this  noise  input  was  too  low  to  appear  on  the  pitch  rate  trace  except 
on  the  most  sensitive  scale,  and  then  is  barely  perceptible. 

If  no  other  frequencies  were  present  in  the  input  to  the  frequency  sensor, 
the  8.5  cps  noise  signal  would  cause  the  gain  to  be  driven  to  the  minimum  per¬ 
mitted  by  the  serve  stop.  By  simulating  pilot  inputs  to  the  elevator,  it  was  pos¬ 
sible  to  work  or  "jack"  the  gain  back  up  to  a  high  value,  but  each  time  the  pilot 
inputs  were  stopped,  the  gain  would  immediately  be  driven  to  the  minimum  value 
by  the  noise  input.  An  elevator  step  input  then  results  in  a  ©  response  closer 
to  the  open-loop  airplane  than  to  the  model. 

Figure  40c  illustrates  the  effect  of  short  bursts  of  noise  on  the  frequen¬ 
cy  s encar  aid  gain  changer.  The  problem  was  started  and  after  about  14  sec¬ 
onds,  a  #tep  c’^vator  input  was  inserted.  The  response  transient  in  the  band- 
passed  signal  indicated  that  the  actuator  mode  had  a  damping  ratio  of  about 
X  =  .  3  -  .4  and  in  the  following  23  seconds,  the  gain  gradually  increased.  At 
this  point,  the  *,  005* /sec,  8.5  cps  noise  signal  was  switched  into  the  rate  gyro 
for  about  4  seconds,  and  the  gain  was  immediately  reduced  toward  the  minimum. 
This  procedure  wat  repeated  two  more  times  on  the  remainder  of  the  record. 

The  recording  of  the  summed  and  filtered  gate  outputs  on  channel  2  of 
Figure  40c  illustrates  the  exponential  rise  of  the  filter  time  constant  when  the 
noise  signal  is  switched  ^rsto  the  circuit.  This  lag  is  also  evident  in  the  re¬ 
sponse  of  the  gain  servo,  ,  ju  channel  3. 

In  order  to  illustrate  the  weighting  effect  of  the  bandpass  filter,  the 
function  generator,  and  the  gate  frequency  :m  noise  inputs,  a  series  of  analog 
records  was  run  with  two  oscillators  feeding  the  frequency  sensor,  one  at 
relatively  high  frequency  and  one  at  relatively  low  frequency.  The  procedure 
used  was  to  set  the  low  frequency  oscillatory  at  a  certain  frequency  (.66  cps  or 
3.5  cps)  with  an  amplitude  of  1. 1  volts  peak  to  peak  and  to  set  the  second  oscil¬ 
lator  at  a  higher  frequency  (8.3,  20  and  35  cps)  anc  adjust  its  amplitude  until 
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the  gain  servo  had  essentially  no  mean  displacement.  This  amounts  to  super¬ 
imposing  the  two  oscillatory  signals  with  the  proper  relative  amplitude  such 
that  the  zero  crossing  time  intervals  of  S9  in  the  frequency  sensor  are  spaced 
so  that  the  gate  output  pulses  due  to  closely  spaced  zero  crossings  integrate  to 
a  value  that  cancels  the  integral  of  the  gate  output  pulses  due  to  zero  crossings 
that  are  far  apart  and  therefore  of  the  opposite  sign. 

Three  sample  records  are  contained  in  Figures  41a,  b,  and  c.  In  Fig¬ 
ure  41a,  the  low  frequency  oscillator  was  set  at  3.5  cps  and  1. 1  volts  double 
amplitude,  and  the  high  frequency  oscillator  was  set  at  8.3  cps  and  1.6  volts 
double  amplitude.  This  combination  of  frequencies  and  amplitudes,  when 
weighted  and  operated  on  by  the  various  components  of  the  frequency  sensor, 
results  in  essentially  no  steady  state  change  in  the  gain  servo  position.  The 
recorder  paper  was  run  at  two  speeds  to  better  show  the  various  frequencies  in 
the  traces.  The  trace  of  the  summed  and  lagged  gate  outputs  is  a  complex  wave 
form  which  is,  however,  fairly  repetitive  or  periodic  with  a  high  frequency  com 
poneni  of  about  2.  3  cps  and  a  low  frequency  component  with  a  period  of  about 
6.3  seconds.  Due  to  *  s  fact  that  the  frequency  of  one  oscillator  signal  was 
about  as  far  above  the  center  of  the  bandpass  filter  as  the  other  was  below  it, 
the  attenuation  of  the  two  signals  was  nearly  the  same  and  the  output  from  the 
bandpass  filter  looks  nearly  the  same  as  the  input  to  it. 

In  Figure  41  b,  however,  the  high  frequency  (20  cps)  input  was  attenu¬ 
ated  much  more  than  the  low  frequency  (3.5  cps)  input  and  the  bandpass  filter 
has  exerted  a  powerful  effect  in  determining  the  zero  crossings  of  the  signal 
$2  .  Again,  the  summed  and  lagged  gate  outputs  appear  to  be  periodic. 

In  Figure  41c,  the  high  frequency  (20  cps)  input  and  the  low  frequency 

(.66  cps), input  are  both  far  from  the  center  frequency  of  the  bandpass  filter 
* 

and  both  are  greatly  attenuated.  Howe/er,  the  very  small  amplitude  signal 
remaining  is  sufficient  to  operate  the  frequency  sensor  ir.  the  absence  of  other 
inputs . 


Figure  4c  is  a  summary  plot  of  the  above  described  data.  In  this  plot. 
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the  ratio  of  the  amplitude  of  the  h'gh  frequency  input  to  the  amplitude  of  the  low 
frequency  input  is  plotted  against  the  frequency  of  the  low  frequency  input  with 
the  frequency  of  the  high  frequency  input  as  a  parameter. 

From  such  a  plot,  it  can  be  determined  whether  the  gain  will  increase 
or  decrease  when  there  are  two  noise  sources  present  of  different  frequency 
and  different  amplitude.  For  example,  if  there  were  two  noise  sources  present, 
one  at  2.0  cps  and  one  at  8.  3  cps,  and  the  amplitude  of  the  high  frequency  input 
was  equal  to  the  amplitude  of  the  low  fre  juency  input,  one  would  expect  the  gain 
to  be  reduced  by  the  frequency  sensor. 

This  study  indicates  that  periodic  noise  in  the  input  signal  to  the  fre¬ 
quency  sensor  can  result  in  the  gain  being  driven  to  one  extreme  or  the  other. 
This  emphasizes  the  importance  of  knowledge  of  structural  modes,  gyro  mount 
location  and  any  other  source  of  periodic  noise  that  may  find  its  way  into  the 
frequency  sensor  input  signal.  The  study  also  suggests  the  possibility  of  using 
other  filters  than  the  simple  bandpass  filter  in  the  frequency  sensor  input  sig¬ 
nal.  It  is  of  inteiest  to  note  that  filters  in  the  frequency  sensor  input  signal 
line  do  not  affect  the  system  root  locus  plots.  Improvement  in  the  performance 
of  the  frequency  sensor  in  the  presence  of  noise  might  also  result  from  a  differ¬ 
ently  shaped  function  generator.  The  present  configuration  results  in  high  fre¬ 
quency  inputs  being  more  effective  in  changing  the  gain  than  low  frequency  in¬ 
puts,  as  illustrated  in  Figure  38.  It  is  suggested  that  this  curve  may  be  re¬ 
shaped  to  advantage. 

Other  possible  ways  to  control  the  effect  of  noise  on  the  frequency  sen¬ 
sor  which  have  occurred  to  the  writer,  but  which  have  not  been  investigated, 
would  be  to  insert  a  very  small  D.  C.  signal  in  to  hold  the  gate  closed  un¬ 
less  the  oscillatory  part  of  St  was  sufficient  to  cause  zero  crossing.  Another 
way  would  be  to  insert  an  oscillator  signal  into  Sg  of  a  frequency  equal  to  (t)&  . 
This  would  dominate  small  amplitude  noise  signals  and  not  result  in  any  gain 
change.  However,  such  an  input  would  have  to  be  very  stable  relative  to  6^ 
and  would  depend  on  precise  gate-open  times  for  it  not  to  affect  the  gain  servo. 
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SECTION  V 
CONCLUSION 


The  intent  of  this  study  was  to  achieve  a  more  thorough  understanding 
cf  the  characteristics  and  behavior  of  self-adaptive  flight  control  systems 
and  the  problemp  that  will  be  encountered  in  designing  such  systems  for  new 
aircraft,  particularly  high  speed  ones.  Three  self-adaptive  concepts,  those 
of  General  Electric,  DODCO,  and  Aeronutronic,  were  evaluated  for  application 
to  re-entry  vehicles  using  the  X-15  airplane  for  the  example.  The  G.E.  sys¬ 
tem  was  studied  in  considerable  depth  and  the  results  are  presented  in  this 
report.  The  DODCO  system,  essentially  a  concept  for  a  nonlinear  gain  element 
employing  switching  logic,  was  also  studied  in  depth  though  its  purely  concep¬ 
tual  status  did  not  permit  the  detailed  type  of  examination  given  the  G.  E. 
system.  The  study  of  the  DODCO  system  is  the  subject  of  Volume  II  of  this 
report  .  The  principal  conclusion  for  this  system  is  that  although  the  nonlinear 
switching  technique  can  produce  effectively  higher  gain  and  hence  smaller 
system  errors  than  in  the  equivalent  linear  system,  the  improvement  is  insuf¬ 
ficient  to  warrant  the  added  complexity  and  difficulties  attendant  with  the  gain 
switching  except  in  very  special  circumstances.  Specifically,  serious  problems 
with  stability  are  encountered  in  high  order  systems,  for  example,  a  flexible 
airframe  with  sensor  dynamics.  The  Aeronutronic  system  was  examined  only 
briefly  because  its  development  was  incomplete  and  significant  examination 
would  have  required  effort  beyond  the  scope  intended  for  this  project.  The 
Aeronutronic  system  concept  involves  the  determination  of  the  system  impulse 
response,  for  which  a  digital  cross-correlator  has  been  successfully  developed, 
and  the  use  of  the  impulse  response  to  adjust  the  system  gains.  The  general 
method,  though  powerful,  lias  a  serious  fundamental  disadvantage  —  a  long 
adaption  time that  is  likely  to  preclude  its  application  to  self-adaptive  flight 
control  systems. 

An  abbreviated  version  of  ASD-TR-61- 104,  Volume  II  can  be  found  in: 

Rynaski,  E. ,  and  Schuler  J. ,  "The  DODCO  Optimal  Control  System  -  A 

Critical  Evaluation",  196?.  Joint  Automatic  Control  Conference,  Paper  12-3. 
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The  conclusions  concerning  the  subject  matter  of  this  report  (Volume  I) 
fall  into  three  categories:  (1)  the  use  of  rnodele  in  high-gain  adaptive  flight 
control  systems,  (2)  the  evaluation  of  the  G.E.  system  as  a  high-gain  self- 
adaptive  system,  and  (3)  the  evaluation  of  the  G.  E.  frequency  sensing  techni¬ 
que  for  gain  adaption.  Though  most  of  the  study  was  concerned  with  control 
of  the  longitudinal  short-period  motions  of  the  airplane  using  pitch-rate 
feedback,  the  results  are  generally  applicable  to  control  of  motion  about  all 
three  airplane  axes  and  a  variety  of  feedback  sensors.  Though  control  system 
and  sensor  dynamics  were  generally  considered,  aircraft  flexibility  was  not. 

High-Gain  Adaptive  Approach 

1.  Reference  models  can  be  used  in  closed-loop  flight  control  systems  to 
tailor  the  response  of  the  aircraft  so  that  it  follows  the  model,  provided 
a  high  enough  system  gain  can  be  realized.  Either  prefilter  models  or 
inverse  models  in  the  feedback  path3  can  be  used  to  achieve  the  desired 
response  to  command  inputs.  Since  disturbances  affect  only  the  closed- 
loop  portion  of  the  system,  the  response  to  disturbances  can  be  tailored 
only  by  inverse  models  in  the  feedback  paths.  Thus  the  response  to 
both  commands  and  disturbances  can  be  tailored,  within  practical  limits, 
by  a  combination  of  prefilter  and  feedback  models. 

2.  Handling  qualities  data  must  be  used  with  care  in  defining  models  for 
control  system  design.  In  the  past  this  data  has  been  misused.  Speci¬ 
fically,  CAL's  longitudinal  data  has  been  inappropriately  applied  to  the 
design  of  pitch-rate  adaptive  systems  by  using  the  normal-acceleration 
model  as  a  pitch-rate  model.  Unsatisfactory  control  systems  have 
resulted.  The  airplane  dynamicist  generating  the  handling  qualities 
data  should  be  more  careful  to  point  out  the  significance  and  the  limi¬ 
tations  of  the  data.  The  control  engineer,  in  turn,  must  become  more 
familiar  with  airplane  dynamics.  It  is  not  enough  just  to  work  from  the 
airplane  transfer  functions. 

3.  Considerable  caution  must  be  exercised  in  the  design  of  high-gain 
systems  to  ensure  that  higher  order  dynamics  (of  order  higher  than 
normally  associated  with  the  uncontrolled  airplane)  introduced  by  th^ 
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control  system  are  properly  taken  into  account.  For  example,  unless 
properly  compensated  for,  an  integrator  in  the  forward  loop,  introduced 
to  eliminate  steady-state  errors,  can  dominate  the  system  response 
or  result  in  an  inappropriate  response  and  so  produce  an  inadequate 
system.  The  effect  of  higher  order  dynamics  can  be  subtle  and  may  not 
be  obvious  in  some  airplane  response  variables,  so  all  important 
response  variables  should  be  examined  to  ensure  that  a  proper  system 
has  been  realized. 

4.  The  high-gain  approach  to  adaptive  flight  control  systems  is  fundamentally 
a  single-input  single-output  approach  (even  though  multiple  loops  are 
often  involved)  and  -tllows  the  use  of  relatively  simple  concepts  of  self¬ 
adaption  to  set  the  system  gains.  In  principle,  the  high-gain  feedback  is 
used  to  eliminate  the  normal  airplane  characteristics  and  desensitize 
the  closed-loop  system  to  changes  in  airplane  parameters;  then  the  desired 
response  is  achieved  by  the  introduction  of  a  model.  This  approach, 
though  feasible,  requires  high  feedback  gains  and  makes  it  difficult  to 
maintain  stability  and  minimize  noise.  The  low-gain  approach,  traditional 
in  stability  augmentation  systems,  employs  only  the  minimum  gain 
required  to  achieve  the  desired  response  since  the  feedback  is  used 
directly  to  modify  the  airplane  parameters.  This  approach  has  the  dis¬ 
advantage  that  the  closed  loop  is  sensitive  to  changes  in  airplane  para¬ 
meters  and  the  gains  must  be  set  precisely,  so  sophisticated  concepts 
of  self-adaption  are  required.  However,  all  the  benefits  realizable  with 
respect  to  stability  and  noise  from  this  use  of  the  low-gain  approach  demand 
that  work  be  continued  on  this  type  of  self-adaptive  system  until  useable 
techniques  are  developed.  The  high-gain  adaptive  systems  should  be 
considered  as  interim  or  sub-optimal  solutions  to  the  problem  of  adaptive 
flight  control  system  design. 

G.  E.  System 

1.  The  G.  E.  system,  which  for  longitudinal  control  employs  a  pitch-rate 

feedback  loop  with  an  inverse  second-order  model, can  produce  a  relatively 
invariant  pitch-rate  response  for  the  longitudinal  short-period  motions 
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for  wide  variations  in  airplane  flight  conditions  and  aerodynamic 
parameters.  However,  the  response  so  produced  is  not  the  response 
desired  for  good  handling  qualities  and  the  pilot  will  not  find  the 
response  invariant.  The  pilot  desires  essentially  an  invariant  second- 
order  response  in  normal  acceleration  and  not  in  pitch  rate.  To 
achieve  the  proper  response  with  the  modeled  pitch-rate  feedback,  a 
prefilter  incorporating  a  variable  gain  and  a  variable  lead  time  constant 

i.e.,  £k  {1  +  s)3  must  be  inserted  between  the  pilot's  control  and 

the  input  to  the  closed-loop  system.  The  variable  gain  (K^)  must  pro¬ 
vide  an  approximately  constant  stick  force  per  normal  accelexation 

(F  / n  ),  so  the  gain  must  vary  inversely  with  airspeed  (V).  The 
8  2 

prefilter  lead  time  constant  (2“^)  must  be  equal  to  the  time  constant 
produced  by  the  integrator  in  the  forward  loop,  and  the  required  pro¬ 
gramming  is  complicated  since  it  is  not  only  a  function  of  flight  condition 
but  also  the  level  of  tne  adaptive  gain.  Because  of  the  difficulties  in 
programming  the  prefilter  for  a  pitch-rate  system,  it  is  recommended 
that  the  use  of  normal  acceleration  (n  )  or  angle  of  attack  (at  )  feedback 
be  considered  for  the  G.  E.  longitudinal  system.  Normal  acceleration 
feedback  would  require  no  programming.  Angle  of  attack  feedback  would 
require  only  gain  programming  with  dynamic  pressure  (or  indicated 
airspeed). 


2.  In  the  present  mechanization,  the  adaptive  gain  is  located  between  the 
forward  loop  integrator  and  the  elevator  actuator.  This  location  causes 
serious  difficulty,  in  the  form  of  spuriouB  transients,  whenever  the 
adaptive  gain  changes  rapidly  and  there  is  a  significant  value  of  trim 
elevator.  It  would  be  advisable  to  mechanize  the  inverse  model  so  that 
the  adaptive  gain  precedes  the  integrator  so  as  to  avoid  the  spurious 
transients  that  are  likely  to  occur  either  in  maneuvering  or  in  low-speed 
flight. 


3.  The  effect  of  applying  pilot  command  inputs  to  the  elevator  actuator 
valve  rather  than  as  pitch- rate  commands  was  found  to  change  the 
basic  response  of  the  system  from  a  pitch-rate  command  system  to  a 
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pitch-attitude  command  system,  with  the  response  amplitude  being 
inversely  proportional  to  the  adaptive  gain.  Such  a  system  would 
probably  be  more  acceptable  to  the  pilot  during  a  re-entry  than  a 
pitch- rate  command  system.  Also,  the  initial  transient  response  of 
the  pitch-attitude  command  system  is  similar  to  the  desired  short- 
period  response.  This  change  is  a  step  in  the  right  direction, so  this 
variant  of  the  G.  E.  system  should  be  examined  to  see  if  the  proper 
system  response  can  be  realized  through  further  modification. 

4.  The  presence  of  a  test  signal  is  fundamental  to  self-adaptive  systems. 
Unless  a  response  is  excited,  no  information  exists  on  which  to  base 
the  adaptive  process.  It  has  been  hypothesized  that  changes  in  charac¬ 
teristics,  such  as  changes  in  transfer  function  parameters  or  poles, 
will  in  themselves  excite  the  adaptive  process.  Careful  reflection 
shows  this  hypothesis  to  be  false.  If  no  explicit  test  signal,  per  se, 

is  generated,  then  either  some  implicit  signal  generator  such  as  a 
limit  cycle  or  noise  is  required,  or  the  pilot  must  continually  generate 
inputs  of  sufficient  magnitude  to  maintain  the  adaptive  process.  The 
lack  of  a  test  signal  in  the  G.  E.  system  appears  to  be  an  undesirable 
feature,  and  consideration  should  be  given  to  providing  some  means  for 
maintaining  the  adaptive  process.  One  possible  method  is  to  provide 
a  small  bias  to  continually  drive  the  gain  up.  Another  one  is  to  provide 
a  monitor  which  would  insert  a  test  signal  whenever  there  had  beeu 
no  appropriate  excitation  of  the  adaptive  loop  for  some  given  period 
of  time. 

5.  Though  the  G.  E.  system  can  handle  changes  in  aircraft  aerodynamic 
parameters  with  ease,  the  system  is  sensitive  to  changes  in  the  control 
system  characteristics.  Basically,  the  frequency  and  damping  of  the 
mode  detected  by  the  frequency  sensor  should  have  a  unique,  well 
defined,  and  accurately  known  relation  to  the  adaptive  gain.  The  gradient 
of  the  detected  mode  frequency  with  gain  must  be  large  for  the  adaptive 
process  to  work  properly.  These  requirements  mean  that  many  of  the 
control  system  characteristics  must  be  known  quite  accurately, 
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especially  those  which  have  their  natural  frequencies  nearest  to  the 
frequency  of  the  detected  mode.  Also,  these  control  system  elements 
must  have  relatively  invariant  characteristics.  In  the  example 
application  to  the  X-15,  the  nonlinear  character  of  the  actuator  caused 
the  adaptive  process  to  misbehave  seriously.  An  instability  was 
encountered  which  resulted  in  a  large-amplitude  elevator  limit  cycle 
with  the  actuator  driving  the  elevator  at  its  maximum  rate,  first  one 
way  and  then  the  other.  As  another  example,  reasonable  values  for 
the  variation  of  the  damping  of  the  rate  gyro  were  assumed,  and  these 
variations  were  found  to  be  intolerable  for  proper  operation  of  the 
adaptive  loop.  The  G.  E.  system  essentially  uses  the  frequency  of 
the  detected  mode  to  set  its  damping  ratio  by  means  of  the  adaptive 
gain,  if  a  damping  ratio  sensor  for  the  detected  .node  were  used 
instead,  then  the  system  would  not  be  as  sensitive  to  changes  and  non- 
linearities  in  the  control  system. 

6.  The  method  of  mechanizing  the  variable  gain  is  important.  In  the  sys¬ 
tem  studied,  the  gain  was  a  linear  potentiometer  driven  by  an  integrating 
servo.  Thus,  a  given  incremental  signal  from  the  adaptive  loop  caused 
a  fixed  incremental  change  in  the  adaptive  gain,  irrespective  of  the 
value  of  that  gain.  This  arrangement  caused  difficulty  because,  when 
the  gain  was  low,  a  small  increment  in  the  adaptive  loop  signal  caused 

a  large  percentage  change  in  the  gain,  resulting  in  a  noisy  system.  The 
adaptive  gain  should  be  a  nonlinear  function  of  the  adaptive  gain  changer 
input  so  that  approximately  equai  percentage  gain  changes  are  realized 
irrespective  of  the  gain  level.  The  nonlinear  function  must  be  selected 
with  care  since  its  form  affects  the  rate  of  adaption. 

7.  The  required  rate  or  speed  of  the  adaptive  process  will  not  generally 
be  determined  by  the  rapidity  of  changes  in  aircraft  flight  condition 
Nonlinearities  in  the  control  system,  nonlinearities  in  airplane  dynamics 
and  changes  in  airplane  parameters  with  rapid  motion  variables  such  as 
angle  of  attack,  and  accommodation  of  the  transients  excited  by  large 
abrupt  inputs  or  noise  will  generally  impose  the  most  severe  conditions 
and  determine  the  required  speed  of  the  adaptive  process. 
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G,  E.  Frequency  Sensor 

1.  The  frequency  sensor  operates  on  zero  crossings  of  the  input  signal 
after  its  various  harmonics  have  been  weighted  by  the  band  pass  filter. 
Its  output,  which  drives  the  gain  changer,  is  weighted  by  the  function 
generator.  In  the  absence  of  a  signal  from  the  detected  mode,  the 
frequency  sensor  will  operate  on  signals  from  lower  and  higher  fre¬ 
quency  modes.  If  the  rate  gyro  picks  up  higher  frequency  modes,  such 
as  structural  modes,  the  gain  can  easily  be  driven  too  low.  Once 

this  occurs  there  is  no  mechanism  in  the  system  for  restoring  the 
gain.  If  low  frequency  modes  are  present,  such  as  low  frequency 
struct"  ••  r* * Om  66  or  autopilot  modes,  then  the  gain  can  be  driven  too 
high.  However,  in  this  case,  the  detected  mode  becomes  unstable 
and  the  resulting  oscillation  will  tend  to  drive  the  gain  back  down  toward 
its  proper  value.  In  the  presence  of  continued  low  frequency  signals, 
the  gain  settles  into  a  nonlinear  oscillation  or  limit  cycle,  accompanied 
by  periodic  bursts  of  the  detected  mode. 

2.  It  has  been  shown  that  improper  operation  of  the  frequency  sensor  and 
gain  changer  is  likely  if  either  noise  is  present  in  the  system  or  the 
control  system,  particularly  the  control  actuator,  is  nonlinear.  The 
improper  operation  can  manifest  itself  as  a  slow  limit  cycle  of  the 
gain,  or  in  severe  cases,  the  gain  can  migrate  to  either  its  upper  or 
lower  stop.  Migration  to  the  lower  stop  could  probably  be  prevented 
by  adding  a  bias  to  the  frequency  sensor  output  which  would  continually 
drive  the  gain  up.  However,  this  bias  would  probably  result  in  a 
continual  limit  cycle  of  the  gain.  The  over-all  operation  of  the  adaptive 
loop,  and  particular  its  sensitivity  to  noise,  could  be  improved  by  the 
addition  of  a  test  signal  to  continually  excite  the  detected  mode.  Mal¬ 
function  of  the  frequency  sensor  due  to  control  system  nonlinearity  can 
be  prevented  by  eliminating  the  nonlinearity  or  its  effect  on  the  detected 
mode.  If  the  effect  of  the  nonline.  -  is  not  eliminated,  then  extreme 
caution  must  be  exercised  in  the  system  design  to  insure  that  the  non¬ 
linearity  will  not  cause  the  gain  changer  to  malfunction.  These 
conclusions  emphasize  the  importance  of  knowledge  of  control  system 
characteristics  and  nonlinearities,  structural  modes,  gyro  mount 
location,  and  any  sources  of  periodic  noise. 
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TABLE  I 

[  FLIGHT  PATH  PARAMETERS  FOR  X-15  RE-ENTRY 

J 


t 

sec 

h 

1000  ft 

M 

mm 

warn 

1 

oc 

deg 

7 

V 

ft/sec 

0 

226 

5.6 

4.7 

C 

0 

5690 

211 

5.3 

5.1 

0 

0 

5690 

■9 

147 

5.5 

9.6 

15 

.4 

5940 

60 

102 

5.9 

26.0 

15 

mm 

5930 

74 

80 

5.4 

32.3 

12.6 

DS 

5200 

90 

77 

4.8 

29.3 

3 

1. 1 

4700 

*  1/2  p  V* 


The  data  were  obtained  from  a  typical  design-altitude  re-entry 
mission  flown  on  the  North  American  Aviation  flight  simulator 
(Reference  8) 


TABLE  II 

NUMERICAL  BATA  FOR  X-15  EQUATIONS  OF  MOTION 


mm 

IS 

** 

l%c 

Ms 

I-# 

■a 

1 JlPf * 

-.0011 

-.0004 

-0. 1569 

-. 1131 

.0016 

.0002 

20 

-.002 

-.0007 

-0.3173 

-0.2317 

.0034 

.0004 

40 

-.0078 

-.0027 

-1.5083 

-1. 1664 

.0114 

.0022 

5-0 

-.0595 

-.0208 

-11. 11 

-8.51 

.112 

.0172 

74 

-. 1541 

-.0539 

-26.41 

-17.5 

.2795 

.0417 

f7 
1  ° 

? 

-. 1322 

-.0463 

-17.  I 

-12.2 

.2767 

.0372 

'T'*tee<s  data  correspond  to  those  of  Table  I. 
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TABLE  III 

LUMPED  PARAMETERS  FOR  &/S&  AND  u/§e 
TRANSFER  FUNCTIONS  FOR  X-15 


t 

sec 

rad/sec 

r 

1/  • 

Ke 

sec*  1 

sec 

■a 

0.396 

IHB m ; 

-.000949 

755 

mm 

0.563 

360 

E9 

1.229 

H 

1 

86.6 

60 

3.33 

11.2 

74 

5. 14 

4.6 

90 

4.13 

.0551 

-.  160 

4.45 

t 

sec 

sec 

0 

-.721 

20 

-.730 

■ 

40 

-.773 

60 

-.766 

74 

-.662 

1 

90 

-.712 

.0031 

These  data  correspond  to  those  of  Table  I  and  result  from 
numerical  substitutions  in  Equations  (2.4),  (2.5),  and  (2.8) 
to  (2.  10). 
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TABLE  IV 

GAIN  (K3)  AND  FREQUENCY  DATA  FROM  ESIAC 
FOR  ACTUATOR  MODE,  <*  =  0  AND  f  =  .  2 


r 

=  .  5  v  sec 

r 

*  ./  ~  sec 

C  8 

O 

=  .2 

■ 

0 

4T  - 

.2 

n 

ESH 

El 

B3F3 

■SB 

■SB 

D 

1700- 

37- 

1100 

30.5 

475 

44 

220- 

36.5 

■ 

1950 

38.5 

230 

20 

1100- 

37- 

620- 

29.7- 

265- 

43.5- 

155- 

I'  l 

1300 

40.5 

730 

31.5 

295 

44.5 

170 

40 

273- 

36- 

170- 

29.5- 

75- 

43.5- 

42- 

36- 

332 

39.5 

193 

31.5 

83 

46 

46 

38 

60 

25.5- 

35.5- 

16- 

29- 

tfl 

43.5 

BSS 

36- 

37.5 

39 

18 

30.5 

Ka 

S  m 

37 

74 

12-  - 

sn 

ESI 

29- 

eh 

42- 

1.7- 

35.5- 

14 

B£9 

Ka 

29.5 

Hfl 

44.5 

1.9 

36 

90 

20 

36.5- 

12 

29 

5.  1 

43.5 

BSR 

35- 

37 

Di 

37.5 
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SYSTEM  AS  INVERSE  MODEL  IN  FEEDBACK 


Fisurt  8  BLOCK  DIAGRAMS  OF  THE  6.E.  SYSTEM 


Figure  9  FUNCTIONAL  DIAORAM  OF  THE  6.E.  FREQUENCY  SENSOR 
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Flgur#  10  COMPLEX-PLANE  PLOT  OF  SHORT-PERIOD  HANDLING 
QUALITIES  DATA  FINN  REFENENCES  IN  AND  IS 
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X-16  AIRPLANE  t  ~  Hue  STEP  COMMAND  INPUT 

UNCONTROLLED  AIRPUNE 


(FORMED  VITN  4  MODEL  AM  NUMERATOR  TINE  CONSTANT) 


r~r  r~t  i  ib  i  i  r~n&~t  m 

TIME  ~  SEC 

Ftfiur*  I!  COMPARISON  OF  9  AMD  &  RESPONSES  FOR  THE  UNCONTROLLED  AIRPLANE 
AND  THE  CONTROLLED  AIRPLANE  WITH  TNO  DIFFERENT  &  MODELS 
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Figure  12  MOT  LOCUS  PLOTS  SHOVINQ  OENERAL  OPERATION  OF  O.E.  SYSTEM 
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CONTROLLED  AIRPLANE  WITH  INPUT-LOOP  PREFILTER  UNCONTROLLED  AIRPLANE 

ST EP  COMMAND  INPUT  STEP  ELEVATOR  INPUT 
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Figure  15  ANALOG  RESPONSES  OF  G.E.  SYSTEM  AND  X-15  AIRPLANE 
AT  VARIOUS  TIMES  DURING  RE-ENTRY  (PART  b) 
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Figure  18  EFFECT  OF  VARIATIONS  IN  ACTUATOR  TIME  CONSTANT 


Figure  19  ATTEMPTED  COMPENSATION  FOR  VARIATION  IN  ACTUATOR  TIME  CONSTANT 
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VARIATION  IN  RATE  GYRO  DAMPING 
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RESPONSE  OF  ADAPTIVE  LOOP  TO  INITIAL  GAIN  OFFSET 
Part  a  Positive  Gain  Offset  -  High  Gain  Case 


Figure  21  RESPONSE  OF  ADAPTIVE  LOOP  TO  INITIAL  GAIN  OFFSET 

fart  a  (continued)  Positive  Gain  Offset  -  High  Gain  Case 


EASE  Analog 
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Figure  2i  RESPONSE  OF  ADAPTIVE  LOOP  TO  INITIAL  GAIN  OFFSET 
Part  b  Positive  Gain  Offset  -  Low  Gain  Case 
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Figure  23  RESPONSE  OF  AIRPLANE  TO  STEP  CHANGES  IN  ADAPTIVE  GAIN 

Part  b  Including  Adaptive- 1  oop  Dynamics 
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K-  =  4  Flight  Condition:  t  -  90 

Actuator  -  Linear  =  0.  1 
Input  -  =  0.3c/sec 

Figure  25  ANALOG  RESPONSES  OF  X- 15  WITH  Q.E.  SYSTEM 
PART  b  RESPONSE  TO  PITCH  RATE  COMMANDS 


ASD-TR-6 1-104,  Vol.  I 


Flight  Condition: 
t  =  0 

Actuator  -  Linear 

7'=  0.1 

Gain:  K  =  150 
Input:  Step  =  26* 


Flight  Condition: 
t=  0 

Actuator  -  Linear 

r  =  0.1 

Gain:  K3  =  250 
Input:  Ste;  §e  =  26* 


REAC  ANALOG 


Figure  25  ANALOG  RESPONSES  OF  X-15  WITH  G.E,  SYSTEM 

FART  c  RESPONSE  TO  STEP  ELEVATOR  COMMANDS,  FLIGHT  CONDITION  t  =  0 
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Figure  25  ANALOG  RESPONSES  OF  X— 15  WITH  G.E.  SYSTEM 

PART  c  (continued)  RESPONSE  TO  STEP  ELEVATOR  COMUNDS 

FLIGHT  CONDITION  t  =  40 


TR-6 1-104,  VoL  I 
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Figure  26d 


ASD-TR-6 1-104,  Vol.  I 


TIME  (SEC) 

DURING  TYPICAL  RE-ENTRY 

Figure  27  GAIN  (X3)  REQUIRED  TO  PRODUCE  £  =  0.02  FOR  THE 
ACTUATOR  POLE  FOR  TWO  ACTUATOR  TIME  CONSTANTS 


«ra> 


ASD-TR-6 1-104,  Vol.  I 


0  02  .04  .06  .08  .10 

Xv  VALVE  TRAVEL-INCHES 


Figure  dO  RATE  OF  CHARGE  OF  FLOW  WiTH  VALVE  DISPLACEMENT 
FOR  NOHUXEAR  ACTUATOR  VALVE 
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ire  36  OPERATION  OF  FREQUENCY  SENSOR  WITH  NONLINEAR  ACTUATOR 
PART  b  FLIGHT  CONDITION  t  =  40 


Figure  37  FUNCTION  GENERATOR  SHAPE  AND  GATE  OUTPUT  FOR  TWO  INPUT  FREQUENCIES 


INFINITE 


FREQUENCY  AND  FUNCTION  GENERATOR  ON 
E  OF  GAIN  SERVO 


K  - 


r: 

t  • 

U 


131 


M  WMeitVtf/  H  4M 


OEG/SEC 


|4«irslai!*)tn«aaamaaaaB«aBa«aBB»aa9aaaaaaaaacaa2aa*| 

it.BSi»*ajtBBBaBBBBBBiBBgBBBBeB)t99»BBBBBBBamBB<l 


aaaaaaaaaaaaeaaaaaaamaaaaaaaaaait 

laaaaaaaaaaaataaaaaaBaaaamaamaStX 


aanaaaaaaaaaaasaaaaaaaaammamamaa 
_  _  C’JHk>B<*aaBBaa>;ft«L.<aaaaaaBaaraac*j« 

WBBBBBBBBfcBBBBBBBB^^P.  -  BBBBBBBBBBBaBB 

‘  BUBBBBBBBBBBB* 

'  ”  riT.V'IfJ'Til *•  .1  %  * = : . •  % ABBBBBB^BBBBBBt 


.ABBBBBB^BBBBl 


OiGAIN  K 


fiaizaaraaaat 


,  , _ ABB  _ 

I  sz?$taaank 

U.'(C«|CU£ 


^aaaaaaraaaaaaaaaaaaaaaaaaa^ataamaaaaaBaSat 


tmaaaammaL s  »*> vaaaaaaaai 


nBmaaaaaaaaaaaaaaaaaaeaaHiKaaaaaaaaaamaal 
ftiSaa£BaMHaBBBBa«aafi»aaa»aaaaBaaaaaagM| 
.  ,CB*»B»eBBBBB«*BB«BBBB»3lr.B  •BCtBBBBBBBBl»«| 

ns>  •>!*■.  't*.f*5*Ba'iBBBi?BBBaaaKk<«eat.*BaBBaaBrBAr 
I**,**.  *<  ^^x--a*',BBaaaai»c;«B^Bv-.vt*safccr«BBB*BaaB 

;  -l *•  .  Of  ..imilBtHtl,;.  >  at-.WBBBB^BBBB? 


.  X 


S. 

VOLTS 


Z*BK 

ta*t 

•  UL 


aaaaatraaaaai 

taaaaaaaaan- 

laaaaaaaaat) 


’rs*>aaaanar 

-« teaamaaaaat 

fl’flfHIBBt 


'»  .«)59i 

rj  •  -1 


MBftt 
t  -S*C^-9SVSi 


iBBaaaaBait 

■IBBCBaSBE* 

nBSBBBBBBB»> 


B&BaaaiaaB 

IBBBOMBBBBI 

m^tHuilV 


IBBlI 


KOtSE"-'  Qh  =  *0.03  OEG/SEC  *  3.5  CPS 


'*K 


pJHH^»aBs.-»B«»__ 

■•aaamama m**r.  ■’•ssbbbbbbb*- 

4anieu*<t«^  iikikicatiMt.  .aii 
*. \i  ra? f/.it .  0Bti44.3*«sBBaBaBB»aa9'i^*BBaBBHBaaaali. 
•!  Vv.  9*  H  -*BC3ilBB<»*C.B**aB»aB*BV»*»B*BBaBiBBBa< 


OEG/SEC 


xaa6iayaa^fBtaaaaae.ss!Jmaaaaaaaaaai* 

t  •***■*  o»*  tmuaatw*3  im&»*aaaaaaaa« 

it-;  tf-jraaum^wa-  mt*a*eaaar  i*i 

i  urirjml.  -xx«i  ««•«*»*<-  ’tick 

.  seSBRlBOW  JIBBlI 

.  p  -  •’  **&*;-* 

.  vi.  uamitii 


ASD-TR-6 1-104,  Vol.  I 


135 


FIGURE  40  EFFECT  OF  NOISE  IN  RATE  GYRO  OUTPUT  ON  OPERATION 
OF  FREQUENCY  SENSOR  AND  GAIN  SERVO 
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cpa  at  1.  1  volts,  <*>,  =  8.3  cps  at  1.6  volta 
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FIGURE  RESPONSE  OF  FREQUENCY  SENSOR  AND  GAIN  SERVO 
TO  INPUT  SIGNAL  COMPOSED  OF  TWO  SINUSOIDS 
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FIGURE  «H  RESPONSE  OF  FREQUENCY  SENSOR  AND  GAIN  SERVO 
TO  INPUT  SIGNAL  COMPOSED  OF  TWO  SINUSOIDS 

Part  c  a]  as  0.66  cps  at  1. 1  volts,  o>t  ~  20  cps  at  2.  1  volts 
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APPENDIX  A 

AIRPLANE  EQUATIONS  AND  PARAMETERS  USED  IN  THE  STUDY 

The  following  equations  of  motion  were  used  to  represent  the  airplane. 
They  assume  that  there  is  no  speed  change,  and  that  the  incremental  effects 
of  gravity  can  be  neglected. 

6  •  +  M^cl  +  MgSe  j 

at  *  Q-L^^-L^Sq 

•  *  . 
jr  x  Q  -  & 

.  The  equations  imply  that  the  reference  axes  are  stability  axes  (wind  axes  if 
no  lateral  motions  are  present)  and  that  the  wings  are  always  level  so  that 

The  reference  condition  for  these  equations  is  not  equilibrium  flight,  but  is 
rather  a  precomputed  re-entry  flight  path.  The  dependent  variables,  (0,0 L  , 
§c,  and  2f  )  are  the  incremental  values  from  the  reference  condition. 

The  following  transfer  functions  arise  from  Equations  (A-l),  and  are  of 
use  in  interpreting  the  analog  solutions  and  root  locus  plots. 


f  (A-l) 

J 


efo)  „ 

5  (Mg  ~L$ Mi.)  ♦  (Mg L*.  * M « f-s) 

(A-2) 

*e(*)  ' 

»*•»  »(*■«-% -«4h-<vv*) 

*(*)  _ 

5  (r^l)  *  (M»  * 

(A-3) 

8e®  ~ 

5^s  Lj 

Ms)  _ 

Ls$*  *  s(-~*  (MsU-M*  LS) 

(A-4) 

SM 

The  natural  frequency  and  damping  in  these  transfer  functions  may  be  expressed 

as: 
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"a  a 

(A-5) 

* 

L*- M+- Mb 

(A- 6) 

Zu)„ 

The  transfer  function!  may  be  written  in  lumped  parameter  form  aa: 


Oils). 

*(*> 


(A-7) 

(A-8) 


TM  B  K*(b*9*+At9+l) 

*•«  ur^w77* 


i* 

-V* 


<A-9) 
(A- 10) 

(A- 11) 


where 


Wjia*  M&Lg 

(A- 12) 

«.* 

Jk*  Pi  Li 

•  Mg  -  M j  4*, 

(A- 13) 

**  * 

jftlzLgMk  ■ 

Mj  L| 

(A- 14) 

**>  * 

Mj  f  tj 

(A- 15) 

*4 

(A- 16) 
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ACn-TO_4  1.IAA 

*  •%-«*“  *  V  *| 


\f^1  T 

f  V4«  * 


A,  = - ** - 

*  Mj  L^^-M^Lg 

a.  -r 


(A- 17) 
(A- 18) 


Flight  conditions  have  been  selected  at  various  times  during  a  typical 
X-15  re-entry  for  the  design  altitude  mission  (Reference  8).  The  flight  condi¬ 
tion  data  for  this  re-entry  are  given  in  Table  I,  the  numerical  data  for  the 
equations  of  motion  in  Table  II,  and  the  lumped  parameters  for  the  and 

the  ***/ Sq  transfer  functions  in  Table  III. 


ASD-TR-6!  ~  104,  Vol.  1 


APPENDIX  B 

ESIAC  COMPUTER  AND  ROOT  LOCUS  COORDINATES 

The  root  locua  technique  is  used  in  this  study  as  the  primary  analyti¬ 
cal  tool  for  handling  linearized  problems.  The  ESIAC  computer  is  used  to 
physically  compute  and  plot  these  root  loci.  This  computer  greatly  facili¬ 
tates  this  work,  but  it  uses  coordinates  which  are  different  from  those  nor¬ 
mally  used.  The  purpose  of  this  appendix  is  to  describe  these  coordinates 
and  show  their  relation  to  the  normal  root  locus  coordinates. 

Description  of  the  ESIAC 

The  ESIAC  solves  the  equation 

Ks"  IE 

1*1  4 

or 

F*  r«"  # 

where  the  5/  are  the  complex  poles  (for^  =  -1,  -2)  and  zeros  (for^  =  +1, 
+2)  of  the  complex  polynomial  F  .  K  is  the  gain. 

The  ESIAC  consists  of  a  console  with  a  plotting  board  at  the  center, 
two  electrical  potential  analogy  sheets  located  on  either  side,  and  controls 
as  shown  on  Figure  43 .  A  frame  positions  moving  probes  on  the  potential 
analogy  sheets. 

The  ESIAC  contains  two  separate  analog  circuits  -  one  for  the  magni¬ 
tude  of  h  ,  If]  ,  and  the  other  for  the  angle  of  P  ,  A.  F  .  In  setting  up  a 
problem,  a  pattern  of  probes  is  placed  in  the  frame  over  each  sheet  at  posi¬ 
tions  corresponding  to  the  's  (the  zeros  and  poles)  in  the  equation.  The 
frames  are  then  moved  over  the  sheets  to  positions  corresponding  to  values 
of  the  complex  variable  s  .  The  voltage  distributions  on  the  sheets  are 
such  that  the  voltage  at  each  probe  represents  on  one  sheet  the  magnitude 
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and  on  the  other  sheet  the  angle  of  one  of  the  (  S  -  s*  )  or  (  1-  */«<•  )  factors. 

The  exponent  JL  on  each  factor  is  established  by  selecting  the  value  of  an 
electric  coupling  element  in  the  clamp  which  attaches  the  probe  to  the  frame. 
Dials  on  the  panel  provide  for  setting  the  values  of  F" ,  td  ,  and  5  ,  and  for 
selecting  the  sign  of  Id  ,  the  value  of  the  exponent  n  ,  and  the  form  of  the 
general  factor  (  S  -  S4')  or  ( 1-  s/s;  ). 

When  the  controls  representing  all  the  parameters  in  the  equation 
are  adjusted  to  any  combination  of  values  ftr  which  the  two  sides  of  the  equa¬ 
tion  are  equal,  a  pair  of  electrical  null  indicators,  one  for  the  magnitude  and 
one  for  the  angle,  show  the  condition  of  equality.  After  the  values  of  all 
given  or  assumed  parameters  have  been  set,  the  equation  is  solved  simply  by 
adjusting  the  value  of  the  unknown  until  null  balances  are  obtained. 

For  root  finding  or  root  locus  plotting,  the  motion  of  the  probe  assem¬ 
bly  is  directly  coupled  to  the  stylus  on  the  plotting  board.  For  automatic 
plotting,  the  probe  assembly  is  moved  to  scan  values  of  5  ,  and  whenever  a 
magnitude  or  angle  balance  is  passed,  the  stylus  automatically  marks  the  po¬ 
sition  on  the  plotting  board.  A  line  joining  the  marks  for  magnitude  balance 
is  the  solution  of  the  magnitude  equation,  and  similarly  a  line  joining  the  marks 
for  angle  balance  is  the  solution  of  the  angle  equation.  The  value  of  5  at  each 
intersection  of  these  lines  is  therefore  a  root  of  the  complex  equation. 

To  use  the  ESIAC  for  frequency  response  plotting,  the  position  of  the 
probe  assembly  in  the  coordinate  representing  d.  S  is  fixed  at  90*,  and  its 
motion  in  the  Is  I  direction  is  coupled  to  the  plotting  assembly  in  the  hori¬ 
zontal  direction.  The  vertical  motion  of  the  plotting  stylus  is  coupled  to  the 
IFI  dial  for  gain  plotting  or  to  the  4.  F  dial  for  phase  plotting.  When  the  cir¬ 
cuit  is  balanced  at  various  values  of  Is  I  ,  the  stylus  plots  the  log  gain  or 
phase  angle  versus  log  frequency. 

Since  the  plotting  stylus  is  directly  coupled  to  the  dials,  it  is  not 
necessary  to  read  the  dials  while  making  these  plots.  The  stylus  is  simply 
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I 

I 

moved  about  over  thi  plotting  paper  and  the  automatic  marking  circuit  plots 
the  curves. 

Coordinates  for  Root  Locus  Plots 

The  polar  coordinates  are  the  magnitude,  / S /  ,  and  angle,  Z  S  ,  of 
the  complex  variable  5  .  The  ESIAC  coordinates  for  the  root  locus  plots 
are  a  transformation  of  the  polar  coordinates  of  the  S  *  d+jai  plane.  The 
transformation  is  shown  in  the  accompanying  three -step  diagram.  Figure  44  . 
The  horizontal  scale,  00 n  ,  is  logarithmic.  The  vertical  scale  is  in  degrees 
with  corresponding  damping  ratios  at  cos ~  1-0. 1,  cos'1-0.2,  ...  etc. 

On  Figure  45  are  root  locus  plots  of  the  function, 

,  ■■ 

V~O.Z2+O.ZOj)0~  022- 0.20 j  ) 

in  both  polar  and  ESIAC  coordinates.  Note  that  the  pole  at  the  origin  is  not 
shown  on  ESIAC  coordinates;  however,  the  information  is  programmed  in  the 
computer  electronically. 
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FIGURE  43  THE  FLIGHT  RESEARCH  DEPARTMENT  ESIAC 
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APPENDIX  C 

EASE  ANALOG  COMPUTER  SCHEMATIC  DIAGRAMS 


ANALOG  SCHEMATIC  -  AIRCRAFT 


FIGURE  47  ANALOG  SCHEMATIC  -  CONTROL  SYSTEM 


+  *0 


Figure  48  ANALOG  SCHEMATIC  -  FREQUENCY  SENSOR 
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TABLE  VI 


COEFFICIENT  POTENTIOMETER  SETTINGS  FOR 
ANALOG  SIMULATION  OF  X-15 

These  potentiometer  settings  are  for  the  analog  simulation 
in  Figure  10,  and  correspond  to  the  data  of  Table  I. 


shown 


* 

u* 

=  0 

L* 

t  =  20 

i 

U  L 

t  =  40 

U  L 

t  =  60 

U  L 

t  =  74 

U  L 

t  =  90 

U  L 

Load 

P50 

.0055 

.0055 

.010  .010 

.039  .039 

.298  .299 

.771  .774 

.661  .664 

5M 

P51 

.005 

.005 

.009  .009 

.034  .035 

.260  .280 

.674  .727 

.579  .633 

.  25M 

P52 

.0785 

.0785 

. 159  .  159 

.754  .757 

0 

0 

0 

5M 

P53 

.  113 

.  126 

.232  .278 

0 

0 

0 

0 

.  1M 

P54 

0 

- 

0 

0 

.278  .301 

.660  .713 

.428  .470 

.  25M 

P55 

0 

- 

0 

.058  .062 

.426  .532 

.875  .932 

.610  .730 

.  1M 

P56 

.  160 

.  160 

.  340  . 340 

0 

0 

0 

0 

10M 

P57 

.020 

.020 

.040  .040 

.220  .220 

0 

0 

0 

10M 

P58 

0 

- 

0 

.036  .037 

.280  .303 

.699  .750 

.692  .744 

.  25M 

P59 

0 

- 

0 

0 

. 172  .  174 

.417  .427 

.372  .381 

1M 

*U  =  Unloaded,  *L  =  Loaded 

NUMERICAL  DATA  FOR  FUNCTION  GENERATOR 
IN  THE  FREQUENCY  SENSOR 


- 37.5  V 


INPUT  VOLTS 
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appendix  d 

REAC  analog  computer  schematic  diagrams 
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FIGURE  50  ANALOG  SCHEMATIC  OF  AIRCRAFT  EQUATIONS 


K  AND  RATE  GYRO 


